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A.  Description  of  the  ASAP  Data  Assimilation  Procedure 


The  AFGL  Statistical  Analysis  Program  (ASAP)  Data  Analysis  Procedure  (DAP) 
was  designed  to  produce  sequential  forecasts  and  analyses  of  the  global 
atmosphere.  Originally,  the  ASAP  analysis,  patterned  after  the  NMC  optimum 

interpolation  (01)  analysis  detailed  by  Bergman  (1979)^,  McPherson  et  al. 

2  3  4 

(1979)  ,  Kistler  and  Parrish  (1983)  ,  and  Dey  and  Morone  (1983)  ,  was 

developed  to  provide  initial  conditions  for  the  AFGL  global  spectral  model 

(GSM)  described  by  Brenner  et  al.  (1982).^  The  analysis  procedure  was 

designed  to  work  with  the  GSM  and  the  Machenhauer  non-linear  normal  mode 

initialization  (NMI)  described  by  Ballish  (1980)^  to  provide  accurate  model 

forecasts. 

In  order  to  produce  a  forecast  for  a  desired  date,  an  analysis  must  be 
available  for  the  starting  date  of  the  model  integration.  This  analysis  is  a 
blend  of  observations  with  a  previous  forecast,  which  in  turn  depends  on  a 
previous  analysis.  Thus,  before  generating  any  long-term  forecasts  for  a 
particular  set  of  dates,  sequential  analyses  based  on  short-term  forecasts 
must  be  available  as  initial  conditions.  These  analyses  are  produced  by  a 
data  assimilation  procedure  consisting  of  a  short-term  (say,  6-hour)  forecast 
and  an  analysis,  followed  by  an  initialization,  another  forecast,  and  so  on. 

The  overall  structure  of  the  ASAP  DAP  is  illustrated  in  Fig.  1.  Table  1 
lists  the  DAP  computer  codes  by  their  permanent  file  name  (PFN)  and  indicates 
which  computer  they  are  presently  designed  to  operate  on. 

1.  Bergman,  K.  H.,  1979:  Multivariate  analysis  of  temperature  and  winds  using 
optimum  interpolation.  Mon .  Wea .  Rev . .  107 .  1423-1444. 

2.  McPherson,  R.  D.,  K.  H.  Bergman,  R.  E.  Kistler,  G.  B.  Rasch,  and  D.  S. 
Gordon,  1979:  The  NMC  operational  global  data  assimilation  system. 

Mon. Wea.  Rev..  107,  1445-1461. 

3.  Kistler,  R.  E.,  and  D.  F.  Parrish,  1983:  Evolution  of  the  NMC  data 
assimilation  system:  September  1978-January  1982.  Mon .  Wea .  Rev . .  110. 
1335-1346. 

4.  Dey,  C.  H. ,  and  L.  L.  Morone,  1983:  Evolution  and  performance  of  the 
National  Meteorological  Center  Global  Data  Assimilation  System: 
January-December  1982.  Submitted  to  Mon .  Wea .  Rev . 

5.  Brenner,  S. ,  C.  H.  Tang,  and  S.  Y.  K.  Yee,  1982:  The  AFGL  Spectral  Model  of 
the  Moist  Global  Atmosphere:  Documentation  of  the  Baseline  Version.  AFGL- 
TR-82-0393,  Air  Force  Geophysics  Laboratory,  Hanscom  AFB,  MA,  ADA129283. 

6.  Ballish,  B.  A.,  1980:  Initialization.  Theory,  and  Application  to  the  NMC 
Spectral  Model.  Ph.D.  Thesis,  UniverTity  of  Maryland,  151  pp. 
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Table  1.  DAP  Computer  Codes 


PFM 

Computer 

Function 

GUM 

Cray 

Machenhauer  Nonlinear  Normal 
Node  Initialization 

FF 

Cray 

Global  Spectral  Model 

PRESPOI 

PRESFPRSOI 

Cray 

AFGL  Cyber 

Converts  output  from  FF  to 
grid  point  form 

SPX 

SFPRSOI 

Cray 

AFGL  Cyber 

Performs  the  surface 
pressure  01  analysis 

PREASAP 

PREASAPOI 

Cray 

AFGL  Cyber 

Vertically  interpolates 
grid  point  T,  u,  v,  q  to  new 
sigma  layer  pressures  and 
converts  these  to  spectral 
form,  then  rederives  grid 
point  values 

ASAPOI 

ASAPOI 

Cray 

AFGL  Cyber 

Performs  the  upper  air  01 
analysis  (including 
calculation  of  residuals  and 
quality  checking) 

POSTASAPOI 

POSTASAPOI 

Cray 

AFGL  Cyber 

Converts  analyzed  grid 
point  values  to  spectral  form 

B.  Preparation  of  Inputs  for  DAP  Components 

The  first  step  in  preparing  for  an  assimilation  sequence  is  preparing  the 
files  that  contain  the  observations.  Both  the  surface  pressure  analysis  and 
the  upper  air  analysis  require  as  input  a  file  with  all  observation  types  on 
it.  A  description  of  the  file  created  from  extracting  desired  data  from  the 
FGGE  II-B  observation  tapes  is  given  in  Appendix  A(2).  The  program  used  to 
create  this  file  on  a  tape  is  BFRDLS/UN=1}0RQUIS,  and  is  used  and  stored  on  the 
AFGL  Cyber.  This  program  is  described  in  detail  in  Appendix  A(l) .  At 
present,  the  code  is  set  up  to  create  a  tape  of  extracted  FGGE  II-B  data  for 
one  date  and  time  only.  When  a  set  of  tapes  containing  the  extracted 
observations  for  all  analysis  times  in  the  assimilation  sequence  is 
constructed,  it  is  mailed  to  the  following  address  for  storage  in  the  Air 
Force  Weapons  Laboratory  (AFWL)  computer  tape  library: 
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Computer  Tape  Library 

AFWL/ADOC 

Bldg.  402 

Kirtland  AFB,  NM  87117. 

If  the  analysis  segment  of  each  cycle  is  to  be  run  on  the  AFGL  Cyber,  the 

extracted  observations  (output  from  BUFF)  can  be  stored  on  permanent  file 

instead  of  tape  if  desired.  When  the  tapes  have  been  received  by  AFWL,  the 

following  job  stream  can  be  executed  on  the  AFWL  Cyber  to  create  permanent 

files  on  the  AFWL  Cyber  which  can  be  “acquired”  by  the  Cray  for  use  in  the 
analyses: 

j  obname , STMFB , MTl , T500 . 

ACCOUNT, username, charge  code-HAN, AFGL, users  phone  #. 

REQUEST, TAPEl ,S, PE, VSN=tape  name  or  number. 

REQUEST , TAPE2 , *PF . 

BEGIN , STOS I , TACLIB , RB=6  4 . NOFILE  =10 . 

CATALOG, TAPE2,pfn,I0=user  identification. 

All  COC  6600  commands  are  upper  case,  whereas  all  Cray  commands  are  lower 
case.  In  this  job  stream,  the  charge  code  is  used  to  login  to  the  AFWL  CDC 
6600:  LOGIN,  (login  sequence  identification).  The  AFGL  tape  name  or  number 
is  the  name  or  number  affixed  to  the  tape  to  identify  it  when  sent  to  AFWL. 
Procedure  STOSI  converts  "stranger”  tapes  to  SCOPE  Internal  files  for  use  in 
the  AFWL  COC  system.  Note  that  since  the  output  tape  from  BUFF  contains  10 
logical  files,  NOFILE  is  set  equal  to  10  for  this  application. 

Once  these  files  are  created  on  the  AFWL  COC  6600,  each  one  may  be 
transferred  to  the  Cray  by  means  of  the  "acquire”  command  when  logged  into  the 
Cray:  acquire  dn=LFN,pdn=PFN,mf=fa,id=USER  I0ENTIFICATI0N,df=cb,  where  LFN  is 
the  logical  file  name  you  wish  to  refer  to  the  file  by  on  the  Cray,  PFN  is  the 
permanent  file  name  on  the  COC  6600,  and  the  USER  lOENTIFICATION  is  the  same 
as  the  one  you  cataloged  the  file  under  on  the  COC  6600.  You  may  then  store 
each  file  on  the  Cray's  mass  storage  system  (like  cataloging  it)  once  it  is 
"acquired”  by  using  the  command:  mass  save  LFN:KASSNAME  where  LFN  is  same  as 
above,  and  MASSNAME  (both  typed  in  lower  case)  is  the  name  of  the  file  as  it 
is  stored  in  "mass". 

If  a  new  data  assimilation  sequence  is  desired  for  a  date  and  time  for 
which  no  previous  sequence  has  generated  an  analysis,  a  starting  analysis  is 
required.  In  the  recent  ASAP  Oata  Assimilation  Experiment  (Halberstam  et  al.. 
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1984)^,  a  FGGE  III-A  analysis  was  used  as  a  starting  point.  A  series  of 
five  AFGL  Cyber  computer  codes,  detailed  in  Appendix  B,  is  used  to  extract  the 
144  X  73  (2.5*  latitude  longitude  grid)  grid  point  set  of  data  at  15  mandatory 
levels  for  Z,  T,  u,  v,  and  RH  and  convert  them  to  30  wave  rhomboidal  (31  x  31) 
spectral  coefficients  of  T,  n  (absolute  vortlcity),  D  (divergence),  and  q 
(specific  humidity)  on  the  12  sigma  layers,  and  of  model  terrain  surface 
pressure.  The  spectral  truncation  as  well  as  the  number  and  distribution  of 
the  model  sigma  layers  can  be  adjusted  easily  in  the  codes  to  match  the  model 
resolution.  The  resulting  file  of  spectral  coefficients,  written  in  binary, 
must  be  converted  from  the  form  of  the  output  of  program  GPTOSPTUVQPS  (the 
last  in  the  series)  to  the  form  required  by  the  initialization  code  GNM.  This 
is  accomplished  through  the  use  of  program  ASTOGS,  included  as  the  sixth  code 
described  in  Appendix  B.  The  30  wave  spectral  coefficients  of  the  terrain 
height  field  are  added  on  as  the  last  record  of  the  output  from  ASTOGS.  This 
output  file,  in  coded  form,  can  be  sent  to  AFWL  on  tape  (in  which  case  it 
should  be  written  onto  the  tape  in  coded  form,  using  a  FORMAT  (IS,F7 .2,A4,A8) 
for  the  first  WRITE  statement  in  ASTOGS,  and  FORMAT(4(1PE20. 13) )  for  the 
others;  then,  on  the  AFWL  Cyber  it  can  be  transferred  to  a  permanent  file 
using  the  COPYCF  utility) ,  or  it  can  be  transmitted  to  the  AFWL  Cyber  in 
binary  form  using  the  NILNET  facility.  If  the  latter  procedure  is  used,  the 
transmitted  file  must  be  converted  to  coded  form  before  being  "acquired"  on 
the  Cray  for  use  in  GNM.  This  can  be  done  using  a  simple  program  on  the  AFWL 
Cyber  in  which  each  record  is  read  in  binary  and  written  using  the  appropriate 
formats  mentioned  above.  The  coded  output  file  is  then  ready  to  be 
transmitted  to  the  Cray  using  the  "acquire"  command,  and  can  be  used  directly 
as  input  to  GNM  to  start  the  data  assimilation  sequence. 


7.  Halberstam,  I.  M. ,  C.  Johnson,  D.  C.  Norquist,  and  S.-L.  Tung,  1984:  Two 
Methods  of  Global  Data  Assimilation  .  AFGL-TR-84-0260,  Contract  F19628- 

82-C-0023,  Systems  and  Applied  Sciences  Corporation. 
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C.  Running  the  Data  Assimilation  Procedure  on  the  Cray 

Procedures  for  running  a  complete  cycle  (CJREP)  and  a  series  of  cycles 
(DNQOI)  have  been  developed  and  are  listed  and  discussed  in  Appendix  C. 

Briefly,  procedure  DNQOI  initiates  the  procedure  CJREP  repeatedly  using 
sequential  dates  and  times  ddtt  where  dd  is  the  date  and  tt  is  the  Greenwich 
hour  of  day.  The  code  ddtt  for  both  initial  and  final  times  of  the  forecast 
in  a  given  cycle  are  indicated.  This  lets  CJREP  know  which  cycle  it  is  to 
perform.  The  observation  file  for  the  analysis  as  well  as  the  form  of  SPX  and 
ASAPOI  for  the  designated  date  and  time  (both  of  these  codes  have  date  and 
time  checks  within  them,  so  date  and  time  must  be  set  within  the  codes,  and 
each  of  these  stored  separately  in  "mass**)  must  have  as  the  last  four 
characters  of  its  name  in  **mass**  the  proper  ddtt  values.  The  output  files 
generated  with  the  cycle  (INIT<1>,GSM<2>,  etc.)  are  then  assigned  these  four 
characters  where  <1>  refers  to  the  first  ddtt  in  each  command  of  DNQOI,  while 
<2>  refers  to  the  second  ddtt  code.  The  DNQOI  procedure  will  need  to  be 
changed  for  each  new  data  assimilation  sequence  (and  versions  of  SPX<2>  and 
ASAP<2>  will  have  to  be  created  and  stored  for  each  sequence),  but  CJREP  will 
remain  the  same  for  all  cases  (unless  the  last  sequence  stopped  before 
completing  a  cycle).  The  DNQOI  procedure  is  put  into  action  using  the 
'*submit**  command  (assuming  the  desired  form  of  DN(}OI  resides  on  local  file 
dnqoi):  submit  i=dnqoi,t=MH,cl=A,pr=S  where  HN  is  the  number  of  minutes 
required  to  run  all  of  the  cycles  in  DNQOI  (the  sum  of  the  number  of  minutes 
required  for  each  cycle,  indicated  by  the  last  number  in  each  "cosmos"  command 
in  DNQOI)  and  A  and  S  indicate  the  desired  class  and  priority  for  the  job 
respectively  (see  AFWL  "BATCH”  manual). 

A  description  of  the  seven  codes  that  make  up  the  data  assimilation 
procedure  (see  Table  1)  is  given  in  detail  in  Appendices  D  and  E.  Each  of  the 
codes  as  it  is  used  in  CJREP  has  several  inputs  and  at  least  one  output.  At 
present,  CJREP  is  constructed  in  such  a  way  as  to  preserve  six  output  files 
from  each  cycle: 

init<l>  initialized  fields  for  the  ddtt  of  the  forecast  initial 

conditions 

gsm<2>  forecast  fields  for  the  ddtt  of  the  forecast  time 

mmmcor<2>  corrected  grid  point  values  of  T,  u,  v,  q  for  ddtt  of  forecast 

analyzed  corrections  and  estimated  analysis  error  for  Z,  u, 

V,  q  at  grid  points  for  ddtt  of  forecast 


coreae<2> 


ninin<2>  analyzed  fields  for  the  ddtt  of  the  forecast 

postout<2>  spectral  root  mean  square  values  for  analyzed  fields  for 
ddtt  of  forecast 

Several  other  files  are  stored  in  "mass"  during  the  course  of  each  cycle,  but 
since  they  are  not  stored  with  a  ddtt  indicator,  they  are  overwritten  in  the 
following  cycle.  The  six  files  listed  above  were  retained  because  of  their 
importance  in  assessing  the  quality  of  the  data  assimilation  procedure.  A 
procedure  (CJXXX)  was  developed  to  store  four  of  the  files  (init<l>,  gsm<2>, 
coreae<2>,  and  nmin<2>)  for  each  cycle  on  a  tape  in  coded  form.  Files  init<l> 
and  coreae<2>  are  in  binary  form,  so  CJXXX  calls  and  runs  programs  (nmrw  and 
mmx  respectively)  which  re-write  the  files  in  coded  form.  Then  CJXXX  gets  and 
executes  sis  and  sisin .  programs  written  by  an  AFWL  consultant,  which  are 
available  to  all  users,  to  place  the  coded  files  sequentially  on  a  tape.  The 
present  form  of  CJXXX  is  given  below: 

^select  savelog=1844 ,  print  log=log<3> 

*select  printlogshsp 

*interrupt  on  software  error  to  err 

*mass  get  nmrw  tapel : init<l> 

*mass  get  dir=/gcbl669/ctss  sis  sisin 
*cft  i=nmrw,  b=bnrarw 
*ldr  bin=bnmrw,  x=xnmrw 
*xnrarw 

♦destroy  tapel 

♦mass  get  gsm:gsm<2>  tape3 : coreae<2>  nm:nmin<2> 

♦mass  get  nortx 
♦eft  i=mmx,  b=bmmx 
♦Idr  bin=bmmx,  x=xmmx 
♦xmmx 

♦cosmos  i=sisin  with  4  tape5  gsm  tape6  nm  /  5 
♦cosmos  i=sis  with  citos  48  4  dnq<3>  /  5 


Type  4:  Satellite  Sounding  Observations 


These  observations  have  the  form  of  several  records  per  observation.  The 
following  sequence  of  READ  statements  should  be  used  to  access  the  data  on  the 
file  that  contains  Type  4  data: 

READ (1,5) ALAT , ALON , IDA , IHR , IDST , PCCl , PI , PCC2 , P2 , TS , PTROP , TTROP , QCTROP , 
NTYPl , NTYP2 , MTYP3 , NTYP4 , ISCR , IQCI , lOP 
5  FORHAT(1X,2F6.2,3I2,2(F2.0,F3.0) ,F3 . 0 ,F4 . 1 . F3 . 0 ,F2 . 0 ,412 , 211 , 13) 

The  variables  NTYPl ,NTYP2 ,NTYP3 ,NTYP4  indicate  how  many  layers  of  each  type  of 
sounding  level  are  reported.  Since  there  are  two  layers  per  record,  the 
number  of  records  present  for  each  type  is  given  by 

NRECl=(NTYPl+l)/2  there  will  be  NRECl  records  (in  the  form  of  the 
:  1  following  READ  statement)  involving  ITYP=01, 

!  I  followed  by  NREC2  records  (in  the  form  of 

:  !  the  following  READ  statement) 

NREC4  =  NTYP4+l)/2  involving  ITYP=02,  etc. 

The  records  have  the  general  form  (where  ITYP=01,02,03,  or  04): 


READ(1,6)ITYP,QC1,PL1,PU1,V1,QC2,PL2,PU2,V2 
6  FORMAT(1X,I2,2(F2.0,2F5.1,F4.0)) 


The  variables  have  the  following  definitions: 


ALAT:  observation  latitude  to  hundredths  of  a  degree 

ALON:  observation  longitude  to  hundredths  of  a  degree  (in  east  long., 

0*  to  360* 


IDA, IHR:  day  and  closest  hour  of  the  observation 

IDSI:  data  source  index  value  (see  Appendix  A,  Table  I) 

PCCl:  percent  of  effective  cloud  cover  in  first  cloud  layer  (in 

percent) 

PI:  pressure  at  average  cloud  top  of  first  cloud  layer  (mb) 


25 


NOTE:  The  four  digit  numbers  in  parentheses  following  many  of  the 
descriptions  above  are  code  table  numbers  for  tables  in  the  following 
document:  Manual  on  Codes.  World  Meteorological  Organization,  Vol.  1, 
International  Codes,  1974  Edition,  WMO  Publication  No.  306,  available 
at  the  AFGL  Research  Library  desk. 
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PZ: 


IQPZ: 

T: 

IQT: 

DPD: 

IQDPD: 

WD: 

WS: 

IQW: 

ITPPT: 


PT: 

WW: 

CLM 

Cl: 

CBH: 

C2: 

C3: 

IPRT: 

ICPR: 

PR: 

IQPR: 

TCA: 

SCA(I) 


sea  level  pressure,  station  pressure,  or  geopotential  height, 
depending  upon  value  of  KPCI  (pressures  are  to  tenths  of  mb,  but 
geopotential  heights  must  be  multiplied  by  10  to  yield  units  of 
geopotential  meters) 

quality  control  for  pressure/height  value  (see  Appendix  A,  Table 
IX) 

temperature  of  the  observation  to  tenths  of  “C 

quality  control  for  temperature  value  (see  Appendix  A,  Table  IX) 
dew  point  depression  of  the  observation  to  tenths  of  ”C 
quality  control  for  dew  point  depression  value  (see  Appendix  A, 
Table  IX) 

wind  direction  of  the  observation  in  degrees 
wind  speed  of  the  observation  in  meters/ second 
quality  control  for  wind  values  (see  Appendix  A,  Table  IX) 
time  period  indicator  for  precipitation  (indicates  time  period 
over  which  precipitation  value  applies)  (see  Appendix  A,  Table 
XII) 

magnitude  of  precipitation  (mm) 

present  weather  (WW  in  sfc  ob  code)  (4677) 

amount  of  low  or  middle  clouds  (N,  in  sfc  ob  code)  (2700) 

n 

types  of  low  clouds  (C  in  sfc  ob  code)  (0513) 
height  of  cloud  base  ("h”  in  sfc  ob  code)  (1600) 
types  of  middle  layer  clouds  (C  in  sfc  ob  code)  (0515) 

H 

types  of  high  clouds  (C  in  sfc  ob  code)  (0509) 

n 

time  period  indicator  for  pressure  tendency  (see  Appendix  A, 
Table  XII) 

characteristic  of  pressure  tendency  ("a"  in  sfc  ob  code)  (0200) 
magnitude  of  pressure  tendency  (mb) 

quality  control  for  pressure  tendency  (see  Appendix  A,  Table  IX) 

total  cloud  amount  (N  in  sfc  ob  code)  (2700) 

,SCI(I) ,SCH(I) :  significant  cloud  amount  (N  ,  see  2700), 

significant  cloud  type  (C,  see  0500),  significant  cloud  height 

(h  h  ,  see  1677);  there  are  at  most  four  levels  of  cloud 
s  s 

information  in  this  sequence 


I 


L 


L, 


1. 


t 
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Type  3:  Surface  (Land  and  Ocean)  Observations 


These  observations  are  taken  at  surface  stations  (ships,  buoys,  manual 
land  stations,  or  automatic  land  stations).  On  the  file,  the  land  based 
observations  come  first,  followed  by  the  ocean  based  observations.  Because 
extraneous  alphanumeric  characters  are  sometimes  found  in  the  elevation 
variable  field  in  the  ocean  based  observations  (which  should  have  a  zero 
elevation),  two  different  READ  formats  must  be  used  to  access  the  data  on  the 
Type  3  file: 

For  land  based  observations  use 

READ  ( 1 , 4  )  ALAT ,  ALON ,  IDA ,  IHR ,  IDSI ,  KPCI ,  EL ,  PZ ,  IQPZ ,  T ,  IQT ,  DPD ,  IQDPD ,  WD ,  WS ,  I()W , 
ITPPT,PT,TCA,WW,CLM,CBH,Cl,C2,C3,IPRT,ICPR,PR,IQPR,(SCA(I),S(rr(I), 
SCH(I) ,1=1,4) 

4  F0RMAT(1X,2F6.2,3I2,I1,F4,0,F6.1,I1,F5.1,1X,I1,F5.1,1X,I1,2F4.0,1X,I1,I1, 
F4.1,7F2.0,/,1X,I1,I2,F3.1,I1,4(3F2.0)) 

For  ocean  based  observations  use 


READ ( 1 , 4 ) ALAT , ALON , IDA , IHR , IDS I , KPCl , EL , PZ , IQPZ , T , IQT , DPD , IQDPD , WD , WS , IQW , 
ITPPT,PT,TCA,WW,CLM,CBH,C1,C2,C3,IPRT,ICPR,PR,IQPR,(SCA(I) ,SCI(I) , 
SCH(I) ,1=1,4) 

4  FORHAT(2F6.2,3I2,I1,A4,F6.1,I1,F5.1,1X,I1,F5.1,1X,I1,2F4.0,1X,I1,I1,F4.1, 
7F2. 0,11, 12, F3. 1,11, 4(3F2.0)) 


The  variables  have  the  following  definitions: 


ALAT: 
ALON: 
IDA, IHR: 
IDSI: 

KPCI: 


observation  latitude  to  hundredths  of  a  degree 
observation  longitude  to  hundredths  of  a  degree 
day  and  hour  of  the  observation 

data  source  index  value  (31-32  are  land  obs,  33-35  are  ocean 
obs)  (see  Appendix  A,  Table  I  in  FGGE-II  formats  documentation) 
pressure  code  indicator;  indicates  the  nature  of  the  value 
stored  in  PZ  (see  Appendix  A,  Table  X) 
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IDSI: 


ITQC: 


P: 

IQCP: 

H: 

IQCH: 

T: 

IQCT: 

TSIP: 

HSIP: 

NWR: 

lORI: 

ITWl: 

ALTWl: 

ALNWl: 

WOl 

WSl: 

IQCl: 

ITWl, 

IWP: 

ICA: 

HC: 


data  source  index,  specifies  the  device  used  to  secure  the 
observation  (see  Appendix  A,  Table  1) 

indicator  for  type  of  quality  check:  0=no  quality  check, 
l=horizontal  quality  check  made,  2=quality  check  against  certain 
limits  was  made 

pressure  at  flight  altitude  (nearest  whole  mb) 
quality  control  for  pressure  (see  Appendix  A,  Table  IV) 
flight  altitude  (nearest  whole  meter) 

quality  control  on  flight  altitude  (see  Appendix  A,  Table  IV) 
temperature  at  flight  altitude  (nearest  tenth  of  a  *C) 
quality  control  on  temperature  (see  Appendix  A,  Table  IV) 
temperature  at  nearest  standard  isobaric  pressure  surface 
(nearest  tenth  of  a  *C) 

altitude  at  nearest  standard  isobaric  pressure  surface  (meters) 
number  of  wind  observations  Included  in  the  report 
optional  record  Indicator;  O=optional  record  not  included, 

1= included 

type  indicator  for  first  wind  ob  in  record  (0=spot,  l=mean, 
9=unknown) 

latitude  (to  tenths  of  a  degree)  of  the  first  wind  ob  in  record 
longitude  (to  tenths  of  a  degree)  of  the  first  wind  ob  in  record 
direction  (degrees)  of  first  wind  ob  in  record 
speed  (ms  ^)  of  first  wind  ob  in  record 

quality  control  for  first  wind  ob  in  record  (see  Appendix  A, 
Table  IV) 

ALTW2,  ALNW2,  WD2,  WS2,  IQC2:  quantities  corresponding  to  above  for 

second  wind  ob  in  record 

indicator  for  weather  phenomena  observed  (see  Appendix  A,  Table 
VII) 

indicator  for  cloud  amount,  observed  from  base  or  top  (see 
Appendix  A,  Table  VIII) 

height  of  cloud,  observed  from  base  or  top  (meters) 


Type  2:  Aircraft  Observations 


These  observations  are  single  point  observations  taken  by  aircraft  at  a 
known  pressure  altitude  above  sea  level.  The  following  sequence  of  READ 
statements  should  be  used  to  access  the  data  for  each  aircraft  observation  on 
the  file  containing  Type  2  data: 

READ( 1 , 3 ) ALAT , ALOH , IDA , IHR , IDSI , ITQC ,P , IQCP , H , IQCH , T , IQCT , TSIP , HSIP ,  IIWR , 
lORI 

3  FORMAT(1X,2F6.2,3I2,I1.1X.F4.0,I1,F5.0,1X,I1,F4.1,1X,I1,F4.1,F5.0,I2,I1) 

The  variable  NWR  contains  the  number  of  wind  observations  included.  Since 
each  wind  observation  record  (that  follows  the  above  record)  contains  two 
individual  observations,  the  total  number  of  wind  records  is  given  by  the 
FORTRAN  statement  NWREC=(NWR-fl)/2.  They  assume  the  following  form: 

DO  10  N=1.NWREC 

READ  ( 1 , 4  )  ITWl ,  ALTWl ,  ALNWl .  WDl ,  WSl ,  IQCl ,  ITW2 ,  AL1'W2 ,  ALNW2 ,  WD2 ,  WS2 ,  IQC2 

4  FORMATCIX, 2(12, 2F4.1,2F3. 0,11)) 

10  CONTINUE 

If  lORI,  the  optional  record  indicator,  has  a  value  of  zero,  no  record  follows 
the  wind  records,  and  the  end  of  the  report  has  been  reached.  For  IORI=l,  the 
following  record  follows  the  wind  records  as  the  last  record  of  the  report: 

READ(1,5)IWP.ICA,HC 

5  FORMAT  (1X,2I2,F5,0) 

The  variables  have  the  following  definitions: 

ALAT:  latitude  of  the  current  aircraft  position  (to  hundredths  of  a 

degree) 

ALON:  longitude  of  the  current  aircraft  position  (to  hundredths  of  a 

degree) 

IDA, IHR:  date  and  closest  hour  of  the  observation 


WS:  wind  speed  of  the  observation  in  meters /second 

IQCW:  quality  control  for  wind  values  (see  Appendix  A,  Table  IV) 

CC1,CC2  both  variables  have  the  values  99999,  used  to  fill  the  record 

CLOUD  DATA  RECORD  (IDSI=75): 

RBAD(  1 . 3  )  ALAT ,  ALOH .  IDA ,  IHR ,  IDS  I .  CLM ,  Cl ,  CBH ,  C2 ,  C3 

3  FORMAT  (1X,2F6.2,3I2.40X.5F2.0) 

The  variables  have  the  following  definitions: 

ALOT:  observation  latitude  to  hundredths  of  a  degree 

ALOH:  observation  longitude  to  hundredths  of  a  degree 

IDA, IHR:  day  and  closest  hour  of  the  observation 
IDSI:  for  this  type  of  record,  IDSI=75  always 

CLM:  amount  of  low  or  middle  level  clouds  (N.  )  on  code  table  2700 

h 

Cl:  clouds  of  lower  type  (C  )  on  code  table  0513 

La 

CBH:  cloud  base  height  (h)  on  code  table  1600 

C2:  clouds  of  middle  type  (C^)  on  code  table  0515 

C3:  clouds  of  high  tjrpe  (Cj^)  on  code  table  0509 

If  the  user  desires  to  use  both  types  of  data,  use  the  read  statement: 

RBAD< 1 , 4 ) ALAT , ALOH , IDA , IHR , IDSI , P , EL , ITYPE , Z , IQCZ , T , IQCT , DPD , IQCDPD , WD , WS , 
IQCW , CLM , Cl , CBH , C2 , C3 

4  FORMAT(1X,2F6.2,3I2,F5.1,F4.0,I2,F5.0,I2,F4.1,I2,F4.1,I2,2(1X,F3.0), 

I2,5F2.0) 

where  if  IDSl3:ll,12,14,15,16,17,  it  would  be  known  that  CLM, Cl, CBH, C2,  and  C3 
would  be  filled  with  9's,  and 

if  IDSI=75,  it  would  be  known  that  P, EL, ITYPE, Z, IQCZ, T, IQCT, DPD, IQCDPD, 
WD,WS,IQCW  would  be  filled  with  9's. 
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Type  1:  Upper-air  Observations 


These  observations  typically  have  the  form  of  a  sounding.  Thus,  for  each 
observation  there  are  several  levels  (usually  at  particular  pressures)  Where 
the  data  are  observed  and  reported.  The  following  formatted  READ  should  be 
used  to  access  the  data  on  the  Type  1  file: 

UPPER  AIR  DATA  RECORD  (IDSI=:11 .12.14.15 . 16 , 17)  : 

READ ( 1 , 2 ) ALAT , ALOM , IDA , IHR , IDS I , P , FL , ITYPE , Z , IQCZ , T , IQCT , DPD , IQCDPD , WD , WS , 
IQCW,CC1,CC2 

2  FORMATdX,  2F6.2,  312,  F5.1,  F4.0,  12,  F5.0,  12,  F4.1,  12,  F4.1,  12,  2(1X, 
F3.0),  12,  2A5) 

The  variables  have  the  following  definitions: 

ALAT:  observation  latitude  to  hundredths  of  a  degree  (southern 

hemisphere  values  are  negative) 

ALON:  observation  longitude  to  hundredths  of  a  degree  (in  east  long.,  0* 

to  360*) 

IDA,  IHR:  day  and  closest  hour  (24  hour  clock)  of  the  observation 
IDSI:  data  source  index  value,  specifies  the  type  of  sounding 

(rawinsonde,  pilot  balloon,  etc.)  (see  Appendix  A,  Table  I  in  FGGE 

II  formats  documentation) 

P:  pressure  of  the  observation  to  tenths  of  mb 

EL:  elevation  of  the  station  above  sea  level  in  meters 

ITYPE:  type  of  level  (standard,  significant,  etc.)  (see  Appendix  A,  Table 

III  in  FGGE  II  formats  documentation) 

Z:  geopotential  height  of  the  observation  in  geopotential  meters 

IQCZ:  quality  control  for  geopotential  height  value  (see  Appendix  A, 

Table  IV) 

T:  temperature  of  the  observation  to  tenths  of  *C 

IQCT:  quality  control  for  temperature  value  (see  Appendix  A,  Table  IV) 

DPD:  dew  point  depression  of  the  observation  to  tenths  of  *C 

IQCPDP  quality  control  for  dew  point  depression  value  (see  Appendix  A, 
Table  IV) 

WD:  wind  direction  of  the  observation  in  degrees 


Appendix  A(2):  Structure  of  Pile  Containing  Data  Extracted  From 
FGGE  II-B  Data  Tapes  (Output  File  From  BUFF) 

Program  BUFF,  located  on  pf n^BUFFREADLS ,  id=NORQUIST,  cy=2,  creates  a  tape 
containing  10  files  consisting  of  a  header  file  and  nine  data  files.  The  nine 
data  files  contain  the  FGGE  II-B  data  for  the  date  and  time  indicated  in  the 
header  file.  End-of-file  marks  follow  each  logical  file  on  the  tape. 

The  header  file  contains  one  record  which  should  be  read  as  follows: 

READ ( 2 , 1 ) IMO . IDA . ITR , IHOUR 
1  FORMAT  (412). 

This  record  is  followed  by  an  end-of-flle  mark,  which  is  in  turn  followed  by 
one  data  file  for  each  of  the  following  types  of  data  in  this  order:  Type  1, 
Type  2,  Type  3,  Type  4,  Type  5,  Type  6a,  Type  6b,  Type  7,  Type  8.  Each  of 
these  files  is  followed  by  an  end-of-file  mark.  The  various  data  types  are 
described  in  the  document  "Formats  for  the  International  Exchange  of  Level  II 
Data  Sets  During  the  FGGE"  which  is  available  for  loan  from  the  author.  All 
references  to  appendices  in  the  following  data  type  descriptions  are 
references  to  that  document. 


67  Statement  Numbers  JDFI  is  the  index  of  data  type 

(01-08),  so  in  order  to  process 
(output)  a  particular  data  type, 
the  corresponding  statement 
number  in  the  "GO  TO"  should  be 
the  statement  number  of  the 
proper  "UNPACK"  subroutine  call. 
If  a  particular  type  of  data  is 
not  desired,  enter  25  for  its 
statement  label  index  in  the  "GO 
TO". 

OUTPUT  FILES: 

TAPE2  Tape  containing  extracted  FGGB  II-B  data  [this  file  is  explained 
in  detail  in  Appendix  A(2)]. 

WRITTEN  OUTPUT: 

1.  JDSI,  date,  and  times  of  any  data  records  skipped  over  by  value  of 
constant  in  statement  36  of  BUFF. 

2.  DSI  (data  source  index),  latitude,  and  longitude  of  observations 
copied  for  types  1,2,4  (these  can  be  turned  off  in  SUBROUTINE  REPIDR 
by  moving  statement  label  20  to  the  RETURN  statement  and  commenting 
out  the  PRINT  statement). 

3.  Indicators  of  end-of-file  marks  written  to  TAPE2  (note  that  besides 
the  EOF  marks  indicated,  an  EOF  mark  is  written  after  the  date,  time 
header  record;  see  statements  15-17  in  BUFF). 

4.  Latitude,  longitude,  and  estimated  terrain  elevation  of  all  satellite 
sounding  observations  rejected  for  being  over  land  (this  may  be 
turned  off  in  SUBROUTINE  UNPACK4,  statement  37). 

5.  Values  of  IREC,  NRECl,  and  IT  in  SUBROUTINE  UNPACK4  (turn  off  by 
moving  label  99  to  RETURN  statement  and  commenting  out  PRINT* 
statement) . 

(NOTE:  The  values  that  are  written  to  TAPE2  in  each  "UNPACK"  subroutine 

may  be  printed  out  by  removing  comment  indicators  in  column  1  from  the 

appropriate  print  statement.  Caution:  if  this  is  done  for  all  data 

types,  output  will  be  very  thick.) 


IHPUT  FILES: 


TAPBl  Copy  of  FGGE  II-B  data  tape  that  excludes  test  files  and  header 
files  (see  "Formats  for  the  International  Exchange  of  Level  II 
Data  Sets  During  the  FGGE"). 

TAPEIO  144X73  grid  values  of  terrain  height  from  FGGE  Fixed  Field  tape. 

Read  as  one  binary  record  into  TH(144,73),  the  values  are  arranged 
on  73  latitudes  from  90*N  to  90*S  at  2.5**  Intervals,  and  for  each 
latitude  on  144  longitudes  from  0**E  (Greenwich)  to  357. 5‘E  at  2.5* 
intervals.  Storage:  pfn=TBRHT144,  id=NORQUIST,  cy=l  (load  from 
CC2467). 

ITEMS  REOUIRIMG  CHARGE  FOR  DIFFEREMT  TIMES.  RESOLUTlOilS.  ETC.: 

LIME  MUMBER(S)  VARIABLE  OR  COMSTAMT  BXPLAMATIOM 


12 

lYR 

Last  two  digits  of  year  of 
observation  you  want  to  extract 

12 

IMO 

Month  (1-12)  of  observation  you 
want  to  extract 

12 

IDA 

Date  (1-31)  of  observation  you 
want  to  extract 

12 

IHOUR 

Hour  (0-24)  of  observation  you 
want  to  extract  (present  FGGE 
II-B  tapes  have  observations  at 
IHOUR»0,6,12,18  only) 

36 

7300 

This  statement  allows  the  code 

to  skip  over  any  records 
preceding  those  for  the  time 
desired.  The  value  used  to 
replace  7300  can  be  estimated 
from  runs  extracting  data  from 
previous  times  on  the  tape  (see 
"FINAL  VALUE  OF  LOOP=”).  If  no 
such  runs  have  been  done,  set 
this  value  to  1. 


Appendix  Ad):  Program  BUFF 


AIKITKACT :  KCCK  II  H  lu  uued  uu  Input,  and  a  reuLrucLured  tape  in 

unpacked,  formatted  form  is  output.  All  Input  values  read  from  the 
tape  are  unpacked,  but  only  selected  values  are  written  as  output. 
This  selection  is  based  on  what  values  are  used  by  the  surface 
pressure  and  upper  air  analyses  except  in  the  case  of  satellite 
soundings,  in  which  all  observations  over  land  (determined  by 
interpolating  billnearly  terrain  height  from  2.5*  lat.-lon.  grid  to 
observation  location)  are  excluded  from  output  tape,  leaving  only 
ocean  based  observations  for  use  in  the  analysis  procedure. 

AFGL  CYBER  STORAGE: 

BFRDLS/un=NORQUIS  (RECLAIM  from  CC4830) 

AFWL  CRAY  STORAGE; 

None  (code  on  AFGL  Cyber  only) 

NOS  CONTROL  CARDS  FOR  AFGL  CYBER; 

NORQUIS , T=2000 , CM=1 70000 . 

USER , NORQUIS , NORQUIS . 

CHARGE, 1049,6670. 

LABEL , TAPE 1 ,VSN=CC4925,NT, D=GE , F=SI , PO=R , R , L=COPyOFL2  7273. 

GET,TAPE10=TRHT144 . 

DEFINE, TAPE2=F2B0412/CT=PU. 

FTN , SL , R=3 , PL=999999 . 

LGO. 

GET,PRSFPR. 

ROUTE , PRSFPR , DC= IN . 


permanent  files  on  the  tape  (the  limit  is  20  files),  then  catalog  the  output 
as  one  file.  To  get  the  desired  program  off  this  file,  use  the  NOS/BB 
commands  SKIFF  to  skip  over  any  files  before,  then  the  COPyCF  to  a  requested 
permanent  file,  and  then  catalogue  that  file.  It  can  then  be  "acquired"  for 
use  on  the  Cray. 

The  codes  described  in  the  following  appendices  are  described  in  the  form 
used  for  the  data  assimilation  experiment  detailed  by  Halberstam  et  al. 
(1984).  References  to  GBR  in  the  appendices  refer  to  the  technical  report  by 
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Gerlach  (1983),  while  references  to  HJIIT  indicate  the  technical  report  by 
Halberstam  et  al.  (1984). 


8.  Gerlach,  A.  H. ,  ed. ,  1983:  Objective  Analysis  and  Prediction 

Techniques-1983 .  AFGL-TR-83-0333,  Contract  F19628-82-C-0023 ,  Systems  and 

Applied  Sciences  Corporation,  A0A142441. 


In  this  example,  there  are  four  files  being  used  in  sisin  and  sis,  in  which 
tapes  is  the  coded  initialized  field  file,  gsm  is  the  forecast  model  output, 
tape6  is  the  coded  corrections  and  analysis  error  file,  and  nm  is  the  01 
analysis  output.  In  the  second  "cosmos"  line  above,  citos  is  the  copy 
routine,  the  record  block  size  is  set  at  48  (this  should  not  change),  and 
four  files  are  being  copied  to  tape  dnq<3>.  In  each  "cosmos"  run  five 
minutes  is  allowed  in  this  example.  This  CJXXX  procedure  is  called  by  CJTAP, 
a  procedure  designed  to  run  several  cases  of  the  procedure  above 
sequentially.  The  form  of  CJTAP  corresponding  to  the  example  of  CJXXX  given 
above  is: 

^select  savelog=1844 ,  print log=logrun 
^select  printlog=hsp 
*mass  get  cjxxx 

^cosmos  i=cjxxx  with  0712  0800  001  /  7 

*cosmos  i=cjxxx  with  0800  0812  002  /  7 

^cosmos  i=cjxxx  with  0812  0900  003  /  7 


so  that,  as  before  <1>=0712,  <2>=0800,  <3>=001  in  the  first  "cosmos"  for 
example. 

Once  the  tapes  containing  the  desired  output  files  are  written,  they  can 
be  sent  back  by  request  (phone  production  control,  AV  8-244-8195)  and  read  at 
AFGL  in  formatted  form  using  the  following  NOS  control  cards: 

LABEL , TAPEl , VSN=  tapename , NT , 0=PE , F=S , P0=R , LB=KN . 

FILE , TAPEl , BT=C , RT=F , FL=80 . 

It  is  important  to  note  that  the  files  will  appear  on  the  tape  in  exactly  the 
same  form  as  they  were  as  individual  coded  files,  and  with  end-of-file  marks 
after  each  file.  Incidentally,  the  sis  and  sisin  procedures  can  be  used  to 
copy  any  coded  file  onto  a  tape,  including  program  files.  Thus,  it  is  a 
convenient  way  to  create  a  "backup"  tape  of  one’s  programs.  The  procedure 
CJBCK,  in  mass  storage  under  user=1844,  is  an  example  of  storing  several 
program  files  on  tape.  To  retrieve  the  files  from  tape  for  later  use,  run  the 
STOSI  procedure  on  the  AFWL  Cyber  with  NOFILE  set  to  the  proper  number  of 
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PCC2: 


percent  of  effective  cloud  cover  In  second  cloud  layer  (in 
percent) 

P2:  pressure  at  average  cloud  top  of  second  cloud  layer  (mb) 

TS:  surface  temperature  (in  *C) 

PTROP:  tropopause  pressure  (In  mb) 

TTROP:  tropopause  temperature  (in  ‘O 

QCTROP:  index  of  accuracy  for  tropopause  data  (in  mb) 

NTYPl,llTYP2,KTYP3,NTyp4:  number  of  layers  reported  for  each  of  the 
respective  observation  types  (see  values  of  ITYP  below  for 
description  of  types) 

ISCR:  status  of  corrections  to  the  observation  (see  Appendix  A,  Table 

XXXVII) 

IQCI:  quality  control  Indicator  (see  Appendix  A,  Table  XXXVIII) 

lOP:  logical  record  number  of  this  record  (always=2) 

ITYP:  Indicator  of  type  of  sounding  level;  ITYP=01,  layer  thickness 

between  reference  pressure  level  and  standard  pressure  level; 
ITYP=02,  layer  preclpltable  water,  between  a  reference  pressure 
level  and  a  standard  isobarlc  surface;  ITYP^Oa,  layer  mean 
temperature  between  two  non-standard  pressure  surfaces;  ITYP=:04, 
layer  preclpltable  water  between  two  non-standard  pressure 
surfaces 

The  meanings  of  the  remaining  variables  depend  on  the  value  of  ITYP.  They  are 
summarized  in  the  following  table: 


Type  5;  Satellite  Clear  Radiance  Data 

In  the  present  version  of  the  code,  these  data  were  not  used  because  It 
was  not  known  how  they  could  be  used  in  the  global  analysis.  If  a  use  can 
be  determined  for  them,  an  UNPACKS  subroutine  could  be  Included  in  the 
BUFFREAD  program.  An  end-of-flle  mark  has  been  Included  as  a  place  holder 
for  the  Type  5  data. 
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Type  6a:  Satellite  Temperature  and  Wind 

These  observations  are  single  point  observations  taken  by  satellite.  The 
records  include  observation  values  of  temperature  and  wind  speed  and  direction 
derived  from  cloud  motion.  The  following  formatted  READ  should  be  used  to 
access  data  on  the  file  containing  Type  6a  data  values: 

READ  (1.7)  ALAT .  ALON ,  IDA ,  IHR .  IDSI ,  ITSWD ,  P ,  IQCP ,  IQC .  T ,  WD ,  WS ,  IQCWD ,  IQCWS ,  OEVE , 
OEPE , IBHM . PDC . ICC , ISWFP 

7  FORMATdX, 2F4. 1, 412, F3.0, 211, 3F3.0, 2X1, 2F2.0,2X, II, F2. 0,211) 

The  variables  have  the  following  definitions: 


ALAT: 
ALON: 
IDA, IHR: 
IDSI: 

ITSWD: 

P: 

IQCP: 

IQC: 

T: 

WD: 

WS: 

IQCWD: 

IQCWS: 

OEVE: 

OEPE: 


observation  latitude  to  tenths  of  a  degree 
observation  longitude  to  tenths  of  a  degree 
day  and  closest  hour  of  the  observation 

data  source  index  value,  specifies  the  type  of  observation 
(see  Appendix  A,  Table  I) 

type  of  satellite  wind  derivation  (see  Appendix  A,  Table 
XXXTII) 

pressure  (mb)  at  effective  wind  level 

subjective  pressure  confidence  factor  (see  Appendix  A, 

Table  XXX) 

quality  control  indicator  for  satellite  winds  (see  Appendix 

A,  Table  XXXVI) 

temperature  (”C) 

wind  direction  (degrees) 

wind  speed  (ms~^) 

subjective  wind  direction  confidence  factor  (Appendix  A, 
Table  XXX) 

subjective  wind  speed  confidence  factor  (Appendix  A,  Table 
XXX) 

objective  estimate  of  vector  error  (ms~^)  (99  used  for 
unspecified) 

objective  estimate  of  pressure  error  (XIO^  mb)  (99  used 
for  unspecified) 
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Type  6b:  Satellite  Sea  Surface  Temperature 

These  observations  are  single  point  observations  of  sea  surface  (water) 
temperature  taken  by  satellite.  The  following  formatted  READ  should  be  used 
to  access  the  data  on  the  Type  6b  file: 

READ (1,8) ALATl . ALONl , IDA , IHRl , IDSI , SSTl , IQCSSTl , ALAT2 , ALON2 , IHR2 . SST2 , 
IQCSST2 

8  FORMAT(lX,2F4.1,3I2,F4.1,I2,2F4.1,I2,F4.1,I2) 


The  variables  have  the  following  meanings: 


JDFI: 

ALATl: 

ALONl : 

IDA: 

IHRl: 

IDSI: 

SSTl: 
IQCSSTl : 

ALAT2: 

ALON2: 

IHR2: 

SST2: 

IQCSST2 : 


data  type  (here,  JDFI^A) 

latitude  (to  tenths  of  a  degree)  of  the  first  observation  in  the 
record 

longitude  (to  tenths  of  a  degree)  of  the  first  observation  in 
the  record 

date  for  both  observations  in  the  record 

hour  (minutes  are  dropped)*  of  the  first  observation  in  the 

record 

data  source  index  value,  specifies  the  type  of  observation  (see 
Appendix  A,  Table  I) 

sea  surface  temperature  (*C)  of  first  observation  in  the  record 
quality  indicator  of  first  sea  surface  temperature  observation 
in  the  record  (see  Appendix  A,  Table  XXVIII) 
latitude  (to  tenths  of  a  degree)  of  the  second  observation  in 
the  record 

longitude  (to  tenths  of  a  degree)  of  the  second  observation  in 
the  record 

hour  (minutes  are  dropped)*  of  the  second  observation  in  the 
record 

sea  surface  temperature  (*C)  of  the  second  observation  in  the 
record 

quality  indicator  of  second  sea  surface  temperature  observation 
in  the  record  (see  Appendix  A,  Table  XXVIII) 


*Minutes  have  been  dropped  from  Type  6b. 
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Type  7:  Oceanographic  Data  Observations 


These  observations  are  taken  aboard  ships  on  which  oceanographic 
observations  are  being  taken.  The  following  READ  statement  should  be  used  to 
access  data  from  the  Type  7  data  file: 


READ ( 1 . 9 ) ALAT , ALON , IDA , IHR . IDSI , PS , T . TW . WD . WS . SST , ISSTI , PT , SOLRAD 
9  FORMAT  (1X,2F6.2,3I2,F6.1,2F5.1,F3.0,F2.0,F5.1,I2,2F3.0) 


The  variables  have  the  following  definitions: 


ALAT: 

ALON: 

IDA.IHR: 

IDSI: 

PS: 

T: 

TW; 

VTD: 

WS: 

SST: 

ISSTI: 

PT: 

SOLRAD: 


latitude  of  the  observation  to  hundredths  of  a  degree 
longitude  of  the  observation  to  hundredths  of  a  degree 
day  and  hour  (minutes  were  rounded  to  nearest  hour)  of  the 
observation 

data  source  index  value  (see  Appendix  A,  Table  1} 

sea  level  pressure  (mb) 

surface  air  temperature  (*0) 

air  wet-bulb  temperature  (*C) 

surface  wind  direction  (degrees) 

surface  wind  speed  (ms'^) 

sea  surface  (water)  temperature  (*C) 

instrument  used  to  measure  sea  surface  temperature  (see  Appendix 
A,  Table  II) 

precipitation  (in  units  of  0.2  mm)  that  fell  during  6  hour 

period  preceding  synoptic  time 

_2 

solar  radiation  (XIO  langleys/minute) 
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Type  8:  Drifting  Buoy  Observations 

These  observations  are  taken  by  automatic  sensors  and  recorders  on  buoys 
deployed  to  drift  in  the  oceans.  The  following  READ  statement  should  be  used 
to  access  data  on  the  surface  drifting  buoy  (Type  8)  data  file: 

READ  ( 1 , 9  )  ALAT ,  ALON ,  IDA ,  IHR ,  IDS  I .  P ,  IQCP .  T ,  IQCT ,  WD ,  WS ,  IQCW .  IW ,  SST ,  IQCSSiT 
9  FORMAT(1X,2F6.2,3I2,F5.1,I1,F4.1,I1,2F3.0,I1,I2,F4.1,I1) 

The  variables  have  the  following  definitions: 


ALAT: 

ALON: 

IDA.IHR: 

IDIS: 

P: 

IQCP: 

T: 

IQCT: 

WD: 

WS: 

IQCW: 

IW: 

SST: 

IQCSST: 


latitude  of  drifting  buoy  to  hundredths  of  a  degree 
longitude  of  drifting  buoy  to  hundredths  of  a  degree 
day  and  closest  hour  of  the  observation 

data  source  index  value,  specifies  the  type  of  sounding  (see 

Appendix  A,  Table  I) 

sea  level  pressure  to  tenths  of  a  mb 

quality  control  for  sea  level  pressure  (see  Appendix  A,  Table  IX) 
air  temperature  in  *C 

quality  control  for  air  temperature  (see  Appendix  A,  Table  IX) 
wind  direction  in  degrees 
wind  speed  in  meters/second 

quality  control  for  wind  values  (see  Appendix  A,  Table  IX) 
wind  indicator  (i^  in  surface  ob  code,  code  table  1855) 
sea  surface  (water)  temperature  (in  ”0 

quality  control  for  sea  surface  temperature  (see  Appendix  A, 

Table  IX) 
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Appendix  B:  Programs  for  the  Construction  of  a  FGGB  III-A 
Analysis  in  Spectral  Form  on  Model  Sigma  Levels 

PROGRAM  MAME;  LVL3RTV 

ABSTRACT:  Using  analyzed  (FGGE  III-A)  heights  and  temperatures  on  mandatory 

levels  (1000-300  mb)  and  the  144  x  73  grid  (2.5*  lat.-lon.  grid) 
terrain  heights,  calculate  the  terrain  surface  pressure  for  the 
144  X  73  grid  using  the  GETPS  routine  (see  Gerlach,  1983).  Then 
every  other  longitude  value  (144-*72)  is  interpolated  linearly  to 
62  Gaussian  latitudes  (73'>62)  forming  a  72  x  62  grid  of  terrain 
surface  pressure  values. 

AFGL  CYBER  STORAGE: 

RTVPS/un=HORQUIS  (RECLAIM  from  CC4830) 

AFWL  CRAY  STORAGE: 

Hone  (code  on  AFGL  Cyber  only) 

HOS  CONTROL  CARDS  FOR  AFGL  CYBER: 


NORQUIS,CM=155000,T=64. 

USER , HORQUIS , NORQUIS . 

CHARGE, 1049, 66 70. 

FTH,SL,R=3. 

LABEL , TAPE 1 , VSN=CC49 19 , MT , D=GE , F=SI , PO=R , R , L=COPYOFL27 176. 
GETP,TAPESyS. 

LIBRARY, TAPESYS. 

GET,TAPE3=TRHT144. 

DEFIHE,TAPE2=PZ7262. 

LGO. 

GET,RTVTUV. 

ROUTE , RTVTUV , DC=IN . 

EOR 


34 


IMPUT  PILES; 


TAPEI  Copy  of  FGGE  III-A  data  tape  containing  analyses  of  height, 

temperature,  winds,  and  relative  humidity  on  mandatory  levels. 
The  first  three  files  on  the  original  tape  have  been  excluded 
from  this  copy  tape  (see  “Formats  for  the  International  Exchange 
of  Level  III  Data  Sets  During  the  FGGE"). 

TAPE3  Same  as  TAPBIO  for  Program  BUFF  [see  Appendix  A(l)]. 

ITMS  RBQUIRIilG  CHAMGE  FOR  DIFFEREMT  TIMES.  RBSOLUTIOMS,  ETC.: 

LIME  IIUMBER(S)  VARIABLE  OR  COMSTAMT  EXPLAMATIOM 

2-3  62  Number  of  Gaussian  latitudes  as 

2  (GLATS)  an  array  size  must  be  changed  if 

another  resolution  is  desired 


3,35 

72 

33 

KHALF 

36 

MDA 

37 

MHR 

Humber  of  longitudes  as  an  array 
size  and  as  a  value  for  LOU  must 
be  changed  if  another  resolution 
is  desired  (must  divide  evenly 
into  144) 

Number  of  Gaussian  latitudes  in 
each  hemisphere 

Date  of  month  of  desired  analysis 

Hour  of  desired  analysis 
(analyses  are  at  12  hour 
intervals) 


OUTPUT  FILES; 

TAPE2  Contains  following  two  records:  (1)  72  x  62  record  of  terrain 

surface  pressure  (cb),  and  (2)  72  x  62  record  of  terrain  heights 
(m). 


WRITTEN  OUTPUT: 

1.  Record  number,  data  type,  surface  type,  level,  date,  hour,  number  of 
words,  grid  type  of  each  set  of  data  encountered  on  TAPBl. 

2.  Any  parity  errors  or  end-of-file  marks  encountered  on  TAPEI  will  be 
indicated  in  the  output. 
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PROGRAM  NAME:  LVL3RTV 


ABSTRACT:  This  code  extracts  the  72  x  62  grid  values  of  Z,  T,  u,  v  from  the 

144  x  73  grid  on  the  mandatory  levels  (1000-50  mb),  interpolates 
them  to  the  72  longitude,  62  Gaussian  latitude  grid,  then  uses 
PTOSIG  routine  (see  Gerlach,  1983)  to  interpolate  T,  u,  v  to  sigma 
layers . 

AFGL  CYBER  STORAGE: 

RTVTUV/un=NORQUlS  (RECLAIM  from  CC4830) 

AFWL  CRAY  STORAGE: 

None  (code  on  AFGL  Cyber  only) 

NOS  CONTROL  CARDS  FOR  AFGL  CYBER: 


NORQU IS,CM=125000,T=10000. 

USER , NORQUIS , NORQUIS . 

CHARGE, 1049, 66 70. 

FTN,SL,R=3. 

LABEL , TAPEl ,  VSN=CC4919 ,  NT ,  D=GE ,  F=SI ,  PO=R .  R ,  L=COP3rOFL2  7176. 
GETP, TAPES YS. 

LIBRARY, TAPES YS. 

ATTACH,TAPE4=PZ7262. 

DEFINE, TAPE2=TUV7262. 

LGO. 

GET.RTVRH. 

ROUTE , RTVRH , DC=IN . 

EOR 

INPUT  FILES: 


TAPEl  copy  of  FGGE  III-A  data  tape. 

TAPE4  output  file  (see  TAPE2,  above)  from  retrieval  of  72  x  62  grid 
values  of  terrain  surface  pressure. 
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ITEMS  REOUIRING  CHANGE 

FOR  DIFFERENT  TIMES. 

RESOLUTIONS .  ETC . : 

LINE  NUMBER(S) 

VARIABLE  OR  CONSTANT 

EXPLANATION 

2,4,213 

2(GLATS) 

62 

Number  of  Gaussian  latitudes 
in  the  desired  grid  (in  213, +1) 

4,49 

72 

Number  of  evenly  spaced 
longitudes  in  the  desired  grid 
(must  divide  evenly  into  144) 

6,7  (T  only),  10,43 
2,3(PTOSIG) 

12 

Number  of  sigma  layers  desired 
(for  SI  and  POH,  +1;  for  T, 
second  dimension  must  be  >12;  in 
PTOSIG,PO,POH,TBO,UBO,VBO  only) 

48,212,213 

31 

Number  of  Gaussian  latitudes  in 
each  hemisphere 

60-65 

DEL 

Sigma  depth  of  each  sigma  layer 
(l=highest  layer,  KP=lowest  layer) 

50 

MDA 

Day  of  the  month  of  the  desired 
analysis 

51 

MHR 

Hour  of  the  desired  analysis 
(analyses  are  at  12  hour 
intervals) 

OUTPUT  FILES; 

TAPE2  Contains  T,  u,  v  values  on  KP  sigma  levels  for  a  72  x  62  (or 

whatever  the  values  of  LON  and  2*KHALF  are)  grid.  (Ordered  by 
records  of  36  words  each,  for  a  given  latitude  and  sigma  layer,  12 
T's,  12  u’s,  12  v's.) 

WRITTEN  OUTPUT; 

1.  Any  parity  errors  or  end-of-file  marks  encountered  on  TAPEl  will  be 
indicated  in  the  output. 

2.  Record  number,  data  type,  surface  type,  level,  date,  hour,  number  of 
words,  grid  type  of  each  set  of  data  encountered  on  TAPEl. 

3.  Mandatory  level  pressures,  heights,  temperatures,  and  winds  from  two 
selected  points  in  the  analysis,  and  sigma  layer  pressures, 
temperatures,  and  winds  resulting  from  the  vertical  interpolation  of 
these  values. 


PROGRAM  MAKE:  LVL3RTV 


ABSTRACT :  This  code  extracts  the  72  x  62  grid  values  of  T  and  R.H.  from  the 

144  X  73  grid  on  the  six  mandatory  levels  (1000-300  mb)  where  they 
are  available.  They  are  interpolated  to  72  longitudes  on  62 
Gaussian  latitudes.  Values  of  R.H.  for  the  upper  six  mandatory 
layers  are  calculated  using  a  linear  interpolation  in  In  p 
between  R.H.  (300  mb)  and  (50/300)*R.H.  (300  mb)  at  50  mb.  These 
mandatory  level  values  are  then  interpolated  linearly  in  &n  p  to 
sigma  layers. 

AFGL  CYBER  STORATE: 

RTVRH/un=M0RQUIS  (RECLAIM  from  CC4830) 

AFVfL  CRAY  STORAGE: 

None  (code  on  AFGL  Cyber  only) 

NOS  CONTROL  CARDS  FOR  AFGL  CYBER: 


NORQUIS , CM=125000 , T=512 . 

USER , NORQUIS , NORQUIS . 

CHARGE, 1049. 66 70. 

FTN,SL,R=3, 

LABEL , TAPEl , VSN=CC4919 , NT , D=GE , F=SI , PO=R , R , L=C0PY0FL27 1 76 . 
GETP.TAPESYS. 

LIBRARY, TAPESYS. 

ATTACH , TAPE4=PZ7  26  2 . 

DEFINE , TAPE2=RH7 26 2 . 

LGO. 

GET , RTVTUVQ . 

ROUTE , RTVTUVQ , DC=IN . 

EOR 
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IMPUT  FILES: 


TAPEl  Copy  of  FGGE  III-A  data  tape. 

TAPE4  Output  file  from  retrieval  of  72  x  62  grid  values  of  terrain 
surface  pressure. 


ITEMS  REOUIRIMG  CHAHGE  FOR  DIFFEREMT  TIMES.  RESOLUTIONS.  ETC.: 
LIME  MUMBER(S)  VARIABLE  OR  CONSTAMT  EXPLAMATIOM 


2-4  62 

4,48  72 

6  (TS  only),  7,  9,  42  12 

222  (first  one  only), 

223,  233 
2  (PTOSIG) 

47  31 

58-63  DEL 


Humber  of  Gaussian  latitudes  In 
the  desired  grid  (In  GLATS  40 
should  be  62  -  doesn't  affect 
results  In  this  version) 

Number  of  evenly  spaced 
longitudes  In  the  desired  grid 
(must  divide  evenly  Into  144) 

Niunber  of  sigma  layers  desired 
(for  SI  and  POH,  +  1;  for  TS, 
second  dimension  must  be  >12; 

In  PTOSIG,  QBO,  PO,  POH  only) 

Number  of  Gaussian  latitudes  In 
each  hemisphere. 

Sigma  depth  of  each  sigma  layer 
(l=hlghest  layer,  KP=lowest  layer) 


output  FILES; 


TAPE2  Contains  R.H.  values  on  KP  sigma  layers  for  a  72  x  62  (or  whatever 
the  values  of  LON  and  2*KHALF  are)  grid.  (Ordered  by  records  of 
12  words  each,  for  a  given  latitude  and  sigma  layer,  12  R.H.'s.) 


WRITTEN  OUTPUT: 


Same  as  for  LVL3RTV  for  T,  u,  v,  except  having  R.H.  as  output. 


39 


Appendix  E:  The  ASAP  Codes 


PROGRAM  HAMB:  PRESP 

ABSTRACT :  This  code  takes  the  forecast  spectral  coefficients  as  input  and 

(1)  reconfigures  them  in  the  form  of  the  output  from  program 
TBSTLMN,  and  (2)  converts  them  to  grid  point  values  of  T,  u,  v,  q 
on  12  sigma  layers  and  x  62  latitude  longitude 

grid. 

AFGL  CYBER  STORAGE: 

PRSFPR/un=NORQUIS  (RECLAIM  from  CC4830) 

AFWL  CRAY  STORAGE: 


Mass  storage:  prespoi 

On  backup  tape  CJBCK  (AFWL  CYBER):  prespoi  (fourteenth  file) 

HOS  COIITROL  CARDS  FOR  AFGL  CYBER: 


NORQUIS ,T=5I2 ,CM=270000 . 
USER , NORQUIS , NORQUIS . 
CHARGE, 1049, 66 70. 

FTN , SL , R=3 , PL=30000 . 

ATTACH , TAPB1=GSH04 12 . 

GET , TAPE5=GLTWT31 . 

DEFINE , TAPE2=SP04 12/CT=PU . 
DEFINE , TAPE4=GR04 1 2 /CT=PU . 
GET,BIN=BINCFFT. 

HAP, ON. 

LOAD, BIN. 

LGO. 

GET.SFPRAN. 

ROUTE , SFPRAN , DC= IN . 

BOR 


A 


. 
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PROGRAM  MAME:  MAIM 


ABSTRACT:  The  AFGL  global  spectral  model  (see  Brenner  et  al.,  1982^).  The 

code  inputs  the  Initialized  spectral  coefficients  and  generates 
forecast  spectral  coefficients  as  output  for  the  desired  length  of 
forecast . 


AFGL  CYBER  STORAGE: 

pfn=FFC,  id=NORQUIST,  cy=l  (load  from  CC2467;  on  CYBER  for  storage 
purposes  only) 

AFWL  CRAY  STORAGE: 


Mass  storage:  ff  (present  set  up  for  MTIME=30,  a  12  hr  forecast);  on 
backup  tape  CJBCK  (AFWL  CYBER):  ff  (second  file) 

IMPUT  FILES; 

TAPES  Output  from  normal  mode  initialization  (see  output  file  from 
program  GMAIN) . 

TAPE2  See  Tung  (1985),  section  on  program  GNM,  FORTRAN  Unit  16. 

OUTPUT  FILE: 

TAPE?  Forecast  spectral  coefficients  (same  record  structure  as  input 
TAPES) . 
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Appendix  D:  Vormal  Mode  Initialization  and  Global  Spectral  Model 


PROGRAM  MAMB:  GMAIH 

ABSTRACT :  The  Machenhauer  non-linear  normal  mode  initialization  (see 

9 

Ballish,  1980,  and  Gerlach,  1982  ).  The  code  inputs  the 
uninitialized  spectral  coefficients  (as  output  from  program 
CONSPB,  or  for  the  first  cycle,  from  ASTOGS),  and  for  the  desired 
number  of  iterations  and  vertical  modes,  performs  the 
initialization  of  T,  n>  D»  and  These  fields  along  with 

q  and  ac-a  output  in  the  same  format  as  they  were  input. 

AFGL  CYBBR  STORAGE; 

pfn=GNMC,  id=HORQUIST,  cy=l  (load  from  CC2467;  on  CYBBR  for  storage 

purposes  only) 

AFWL  CRAY  STORAGE; 

Mass  storage:  gnm 

On  backup  tape  CJBCK  (AFWL  CYBER):  gnm  (nineth  file) 

IMPUT  FILES: 

TAPES,  TAPE16,  TAPB80  See  Tung  (1985),^^  section  on  program 

GNM,  for  further  details. 

TAPB18  Spectral  coefficients  of  n,  D,  T,  q,  *gfj.  (see 

output  file  for  program  ASTOGS). 

OUTPUT  FILE: 

TAPB19  Initialized  spectral  coefficients  for  input  to  global  spectral 
model  (record  structure  is  identical  to  input  file  TAPB18  (see 
output  file  for  program  ASTOGS)). 


9.  Gerlach,  A.  M.,  ed.,  1982;  Objective  Analysis  and  Prediction 

Techniques ■  AFGL-TR-82-0294 ,  Contract  F19628-82-C-0023,  Systems  and 

Applied  Sciences  Corporation,  AOA131465. 

10.  Tung.  S.-L.,  1985:  Users  Guide  for  Normal  Mode  Analysis  Method  of  Global 
Data  Assimilation.  AFGL-TR-85-XXXX,  Contract  F19628-82-C-0023,  Systems 
and  Applied  Sciences  Corporation  (in  press). 
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53  •xpi  /  3 

34  *coiipact  tap«3  tapt3 
55  aconpact  tapa6  tap«6 

54  «natt  ttora  tap«3:«ptadqp«sl12n  tapt4tap62tl12n8 

57  *Nafs  itori  outputipaout 

58  ^switch  tapt4  tape?? 

5?  *suitch  tape?  tape? 

60  esuitch  tapet  tape88 

61  eswitch  tape2  tape87 

62  *SMitch  tap«3  tapel 

63  ^switch  tape6  tape4 

64  eswitch  tape?8  tape? 

65  esuitch  tape??  tape2 

66  edeftroy  tapeS 

6?  esuitch  tape?6  tapeS 

68  *Nass  set  tape10:estprer  a«:afap<2> 

6?  *cft  i«as,b>afb 

?0  •Idr  bin>(asb,bcfft),xsaix 

?1  *asx  /  25 

?2  ecoNpact  taped  taped 

?3  aconpact  tapel 1  tapel 1 

?4  aconpact  output  output 

?5  aiiass  store  taped:nN(icor<2>  output:asout 

?6  anass  store  tapel 1 :coreae<2> 

77  asuitch  tape4  tape?? 

?8  asuitch  tape??  taped 
??  asuitch  tape2  tape?8 

80  asuitch  tapet  tape2 

81  asuitch  taped  tapel 

82  anass  s^t  postspostasapoi 

83  acft  impost, bspostb 

84  aidr  bin><postb,bcfft) ,xapostx 

85  apostx  /  10 

86  acoepact  output  output 
8?  aconpact  tape6  tape6 

88  anass  store  outputipostout<2>  tape6iRNin<2> 

8?  adestroy  taped 
?0  adestroy  tapel 
?1  adestroy  tapelO 
?2  adestroy  tape?? 

?3  adestroy  tape? 

?4  adestroy  taped 
?5  adestroy  tape88 
?6  adestroy  tape? 

??  adestroy  tape?8 
?8  adestroy  tape5 
??  adestroy  tape? 

100  adestroy  tape8? 

101  adestroy  tapel 1 

102  asoto 

103  a#rr<2>i 


Fig.  C-1.  CJREP  procedure  (Cjnt.) 


1  Mtltct  favtlo9>1844,printla9>lo3<2> 

2  *intcrrupt  on  softMarterror  to  trr<2> 

3  onats  gat  ff  tapcS  taptSO  tape1<&  gnn  tape18:nnin<1> 

4  «cfi  i>9nN,b>snNb 

5  aldr  bin<9nNb,x<>gnNX 

6  agnnx  /  4 

7  aconpact  tapel?  tapa19 

8  atwitch  tape19  tapc8 

9  anass  git  tapi2 

10  anass  store  tapt8:init<1> 

11  acft  i»tf,b«ffb 

12  aldr  bin=ffb,x=ffx 

13  affx  /  10 

14  aconpact  tape?  tape? 

15  anass  store  tape7:gsN<2> 

14  asHitch  tape?  tapel 

1?  adcstroy  tapeS 

18  adestroy  tapel 6 

19  adestroy  tape80 

20  adestroy  tape8 

21  adestroy  tape? 

22  adestroy  taped 

23  adestroy  tape9 

24  anass  get  tapeSiglategtSIc  pstprespoi  cfft 

25  acft  i«cfft,b*bcfft 
24  acft  iBps,b«psb 

2?  aldr  bin>(psb,bcfft),x«psx 

28  apsx  /  5 

29  aconpact  tape?  tape? 

30  aconpact  taped  taped 

31  anass  store  tape2ispted<|pssl1?8  tape4:spgp4?sl128  output:psout 

32  aiMitch  tapeS  tape99 

33  adestroy  tapel 

34  asMitch  tape?  tapel 

35  anass  get  sp:apx<2>  tape?iglat42  tape2jfgge2<2> 

34  anass  get  tape4:terrsst31c  tapeSitersstd? 

37  acft  i>sp,baspb 

38  aldr  bina(spb,bcfft),x>spx 

39  aspx  /  IS 

40  aconpact  tape9  tape? 

41  anass  store  tape9isfprscor8  outputispout 

42  anass  get  paipriasap 

43  acft  iapa,babpa 

44  aldr  bina<bpa,bcfft),x>xpa 

45  adestroy  tapel 

44  asuitch  taped  tapel 

47  asuitch  tape99  taped 

48  asuitch  tape?  tapa98 

49  asuitch  tapeS  tape? 

50  asuitch  tape2  tape97 

51  asuitch  tape9  tape? 

52  asuitch  taped  tape94 


Fiy.  C-1. 


C.JKKP  Prof.-o'lurf,- 


Appendix  C:  Procedures  for  Running  Data  Assimilation  Cycles 
on  the  AFWL  Cray  Computer 

I.  CJREP:  A  Procedure  for  Running  a  Complete  Cycle 

Fig.  C-1  is  a  listing  of  the  procedure  CJREP,  in  mass  storage  under 
user=1844  and  also  on  backup  tape  CJBCK  (the  handwritten  additions  have  been 
made  to  the  existing  procedure).  This  procedure  automatically  runs  the 
sequence  of  codes  that  make  up  one  assimilation  cycle:  gnm,  ff,  prespoi,  spx, 
preasap,  asapoi,  postasapoi.  The  commands  used  in  this  procedure  can  be 
thought  of  as  control  cards  for  the  Cray.  The  basic  sequence  that  is  repeated 
for  each  of  the  seven  codes  is:  (1)  the  code  and  all  input  files  are  obtained 
(attached)  from  mass  ("mass  get"),  (2)  the  code  itself  is  compiled,  and  a 
binary  file  of  the  compiled  code  is  formed  ("eft, is"),  (3)  the  binary  file  is 
loaded  into  an  executable  file  (idr,  bin=  ),  (4)  the  executable  file  is  run, 

(5)  all  output  files  are  compacted  (this  is  necessary  because  sometimes  longer 
files  may  exist  in  several  parts),  (6)  output  files  are  stored  on  mass  ("mass 
store"),  (7)  files  needed  by  the  succeeding  program  have  their  logical  file 
name  "switched"  to  the  proper  one  for  that  job,  and  (8)  any  logical  files 
no  longer  needed  are  "destroyed". 

II.  ONQOI:  A  Procedure  for  Running  a  Sequence  of  Cycles 

The  following  listing  of  the  file  01K)0I  is  given  as  an  example  of  how  CJREP 
can  be  called  in  sequence  for  several  contiguous  assimilation  cycles: 

*savelog=1844 ,printlog=logmstr 
*mass  get  cjrep 
*mass  get  cj091 
*cosmos  cj09X  with  0912  /  44 
*cosmos  i=cjrep  with  0912  1000  /  88 

*cosmos  i=cjrep  with  1000  1012  /  88 

*cosmos  i=cjrep  with  1012  1100  /  88 

♦cosmos  i=cjrep  with  1100  1112  /  88 

♦cosmos  i=cjrep  with  1112  1200  /  88 

By  using  the  "submit"  command  (as  explained  in  the  text  of  this  users  guide) 
for  the  file  DNQOI ,  part  of  a  cycle  (on  cj09X)  for  the  9th  at  OOZ  to  the  9th 
at  12Z,  and  five  complete  12h  cycles  following  that  one  can  be  executed.  As 
explained  in  the  text,  the  starting  and  ending  date-time  code  ddtt  indicate 
which  cycle  is  processed,  so  by  adjusting  these  values  and  adding  more  or 
using  fewer  "cosmos"  commands,  the  desired  sequential  cycles  may  be  run. 


IF(NSK.EQ.O)  GO  TO  20 
DO  10  NSKP=1,NSK 
10  REAO(l) 

20  RBAD(l)  ((A(M,N.KV) ,M=1,MS) .N=1.HS) 
IF(NSKPA.EQ.O)  GO  TO  40 
DO  30  NSKP=1.NSKPA 
30  READd) 

40  CONTINUE 
WRITE(3)  A 
50  CONTINUE 

DO  60  NTAPE=1,2 

REA0(NTAPE)((A(M.N,1).M=1.MS).N=1,NS) 

WRITE(3}((A(M.N,1).M=1,MS).N=1.NS} 

60  CONTINUE 
STOP 
END 


OUTPUT  FILE; 

TAPE3  Spectral  coefficients  of  Ot  <1  for  31  waves  (0-30)  and  12 

sigma  layers  in  the  following  form  (each  line  below  represents  a 
record) : 

NSTEP,  TIME,  ITIME,  IDATE 
n(I=l,31;  J=l,31;  K=l,12) 

D(I=1,31;  J=l,31;  K=l,12) 

T(I=1,31;  J=l,31;  K=l,12) 
q(I=l,31;  J=l,31;  K=l,12) 

pg^^(I=l,31;  J=l,31)  (this  is  In  of  surface  pressure) 
^^^^(1=1,31;  J=l,31)  (this  is  geopotential  of  the  surface) 
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PROGRAM  NAME:  ASTOGS 

PURPOSE:  This  code  converts  the  spectral  coefficients  from  the  form  of  the 
output  of  program  TESTLMN  to  the  form  required  as  input  to  the 
non-linear  normal  mode  initialization.  Also,  spectral  coefficients 
of  model  terrain  height  are  added  as  the  last  record  to  the  output 
file. 


HOS  CPHTROL  CARDS  AMD  PROGRAM  LISTIMG: 

NORQUIS . CM=150000 . 

USER . NORQUIS , NORQUIS . 

CHARGE, 1049, 66 70. 

FTN,SL,R=3. 

ATTACH, TAPE1=F06 13 . 

GET , TAPE9=TRSST31 . 

DEFINE , TAPE?  =F3A06 13 / CT=PU . 

LGO. 

PURGE, F06 13. 

EOR 

PROGRAM  ASTOGS ( OUTPUT, TAPE2,TAPE2,TAPE3) 

COMPLEX  A(31,31,12) 

DATA  NSTEP,TIME,ITIME,IDATB/0,0.0,4H  00Z,8H  07JAN79/ 
MS=31 
NS=31 
ML=12 

WRITE(3)  NSTEP,TIME,ITIME,IDATE 
DO  50  MT=1,4 
IFrNT.EQ.l)NSK=l 
IF(NT.EQ.2)NSK=2 
IF(MT.EQ.3)MSK=0 
IF(NT.EQ.4)NSK=3 
NSKPA=3-NSK 
REWIND  1 
DO  40  K=1,ML 
KV=NL-K+1 


TAPB4  31  records  of  the  sine  of  the  Gaussian  latitudes  from  equator  to 

pole,  followed  by  31  records  of  the  corresponding  Gaussian  weights 
(all  values  are  correct  to  20  decimal  places). 


ITEMS  REOUIRIMG  CHAMGE  FOR  DIFFBREMT  TIMES.  RESOLUTIOMS.  ETC.: 


LIME  MUMBER(S)  VARIABLE  OR  COMSTAMT  EXPLAMATIOB 


3,  13,  21 


62 


Number  of  Gaussian  latitudes 


2,  3,  21 

72 

2,  3  (FFTl) 

4,  5,  6,  21 

31 

2,  3  (LEGSUM) 

2,  3  (LMN) 

2  (FFTl) 

2,  4  (PMNS) 

9 

4  (FFTl) 

303 

21 

12 

Number  of  evenly-spaced 
longitudes 

Number  of  waves  In  the  spectral 
rhomboldal  truncation 
(32  =  wave  numbers  0  through  30; 
32  =  one  +  rhomboldal  truncation) 


This  value  must  be  >4*NP-t-15 

Number  of  sigma  layers 


OUTPUT  FILE; 


TAPE3 


31  X  31  spectral  coefficients  of  T,  n*  0,  q, 
following  form  (each  line  below  represents  a 
T(I=1,31;  J=l,31;  K=1 
ti(I=1,31;  J=1,31;  K=1) 

D(I=1,31;  J=l,31;  K=l) 
q(I=l,31;  J=l,31;  K=1 


and  pgfc  In  the 
record) : 


T(I=1,31;  J=l,31;  K=12) 
n(I=l,31;  J=l,31;  K=12) 
D(I=1,31;  J=l,31;  K=12) 
q(I=l,31;  J=l,31;  K=12) 
Pg^^(I=l,31;  J=l,31) 


WRITTEN  OUTPUT:  None 


PROGRAM  MAMK:  TRSTLMII 

I  ABSTRACT :  This  code  converts  gridded  (72  x  62)  fields  of  T,  u.  v,  q  on  12 

sigma  layers  and  converts  thon  to  spectral  coefficients  of  T,  n 
(absolute  vorticity),  D  (divergence)  and  q,  and  does  the  same  for 
natural  log  of  surface  pressure. 

I 

AFGL  CYBER  STORAGE : 

GPTOSP/un=NORQUIS  (RECLAIM  from  CC4830) 

I 

AFWL  CitAY  STORAGE: 


None  (code  on  AFGL  Cyber  only) 


NOS  CONTROL  CARDS  FOR  AFGL  CYBER; 

NORQUIS.CM=120000,T=64. 

USER . NORQUIS . NORQUIS . 

CHARGE, 1049, 66 70. 

FTB,SL,R=3. 

ATTACH, TAPE1=PZ7262. 

ATTACH , TAPB2=TUVQ . 

GET , TAPE4=GLTWT31 . 

DEFINE, TAPE3=F06 13. 

GET,A=BINCFFT. 

LOAD, A. 

LGO. 

PURGE, PZ7262. 

PURGE, TUVQ. 

GET,ASTOGS,DC=IN. 

EOR 

INPUT  FILES: 

TAPEl  72  X  62  field  of  values  of  p  -  from  LVL3RTV  (p  ,  ) . 

*^sfc  *^sfc 

TAPE2  72  X  62  field  of  values  T,  u,  v,  q  for  each  of  12  sigma  layers 
from  RW. 
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r'S 


3,5.29 


OUTPUT  FILE; 


Humber  of  Gaussian  latitudes 
Humber  of  evenly-spaced 
longitudes  (must  divide  evenly 
into  144) 

Humber  of  sigma  layers 


TAPB2  Contains  72  X  62  fields  of  T,  u,  v,  R.H.  on  12  sigma  layers. 

Records  are  ordered:  [T(72,62),  u(72,62),  v(72,62),  q(72,62), 
K=ll , . . . , [T(72,62),  u(72.62),  v(72,62),  q(72,62),  K=12J 

WRITTEH  OUTPUT:  Hone 


PROGRAM  HAME;  RW 


ABSTRACT:  This  code  reads  input  from  code  that  makes  72  X  62  T,  u,  v,  R.H. 

file  and  converts  R.H.  to  specific  humidity  (q)  using  retrieved 
surface  pressure,  forming  a  file  of  72  X  62  values  of  T,  u,  v,  q 
for  12  sigma  layers. 

AFGL  CYBER  STORAGE: 

RHTOQ/un=NORQUIS  (RECLAIM  from  CC4830) 

AFWL  CRAY  STORAGE: 

Hone  (code  on  AFGL  Cyber  only) 

MOS  COMTROL  CARDS  FOR  AFGL  CYBER 

NORQUIS,CM=120000,T=16 . 

USER , NORQUIS , NORQUIS . 

CHARGE, 1049, 66 70. 

FTN,SL,R=3. 

ATTACH , TAPEi=TUVRH . 

ATTACH , TAPE4=PZ7  26  2 . 

DEFINE , TAPE23TUVQ . 

LGO. 

PURGE, TUVRH. 

GET.GPTOSP. 

ROUTE , GPTOSP , DC=IN . 

EOR 

INPUT  FILES: 

TAPEl  Output  from  RW  (T,  u,  v,  R.H.). 

TAPE4  Output  file  from  retrieval  of  72  X  62  grid  values  of  terrain 
surface  pressure. 
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TTEMS  RBOUIRIIIG  CHAMGE  FOR  DIFFEREMT  TIMES.  RESOLUTIOMS.  BTC.: 

LIMB  MUMBBR(S)  VARIABLB  OR  COMSTAMT  BXPLAMATIOM 


2,  20-24 

62 

Number  of  Gaussian  latitudes 

2,  17,  21-24 

72 

Number  of  evenly-spaced 
longitudes  (must  divide  evenly 
into  144) 

7,  28 

12 

Number  of  sigma  layers  (for  line 
28,  +1) 

OUTPUT  FILE: 

TAPB2  Contains  72  x  62  fields  of  T,  u,  v,  R.H.  on  12  sigma  layers. 

Records  are  ordered:  [T(72,  62),  u(72,  62),  v  (72,  62),  R.H.  (72, 

62),  K=ll . [T(72,  62),  u  (72,  62),  v  (72,  62),  R.H.  (72,  62), 

K»12] . 


WRITTBM  OUTPUT:  None 


IMPUT  FILES; 


TAPEl  Output  from  global  spectral  model;  spectral  coefficients  in  format 
identical  to  the  output  from  ASTOGS. 

TAPES  31  Gaussian  latitudes  (in  sine  of  latitude)  from  equator  to  pole 
correct  to  20  decimal  places,  followed  by  31  Gaussian  weights  in 
same  orientation  and  accuracy. 


ITEMS  REQUIRING  CHANGE  FOR  DIFFEREMT  TIMES.  RESOLUTIOMS.  ETC.: 


LINE  NUMBER(S) 

VARIABLE  OR  CONSTANT 

EXPLANATION 

3-5,9,13,15,16,23,24 

2.3  (LEGSUM) 

2(FFT1) 

2.4  (PMNS) 

2-4,8  (SPTOGP) 

2-4,7  (UMVM) 

2(FCST0P) 

31  (and  32  when 
paired  with  31) 

Rhomboidal  truncation  (0-30 
waves)  for  spectral  coefficients 
(32  should  be  changed  the  same 
way  as  31;  in  line  16,  KHALF  is 
number  of  Gaussian  latitudes  in 
each  hemisphere 

6,14,25 

2,3  (REORD) 

3  (FFTl) 

7  (SPTOGP 

62 

Number  of  evenly-spaced 
longitudes  in  the  analysis  grid 

26 

NP 

Number  of  Gaussian  latitudes  in 
the  analysis  grid 

3,8,27 

2,3  (REORD) 

2,7  (SPTOGP) 

2  (UMVM) 

2  (FCSTOP) 

12 

Number  of  sigma  layers 

11 

4  (FFTl) 

263 

Dimension  of  the  array  WSAVB 
(should  be  >4*MP+15) 

OUTPUT  FILES: 

TAPE2  Forecast 

TESTLMN. 

spectral 

coefficients 

in  same  format  as  output  from 

TAPE4  62  X  62  X  12  grid  values  of  T, 

Psfc  arranged  as  follows: 

u.  V,  q,  and  62  x  62  grid  values 

DO  10  LAT=1,62 
DO  10  LON=l,62 

T(LAT,L0N,K=1,12) ,  U(LAT,LOH,K=l,12) ,  V(LAT,L0M,K=1 , 12) 
q(LAT,L0M,K=l,12),  Pgfc(LAT,L0N) 

10  CONTINUE 
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WRITTEM  OUTPUT; 


1.  For  the  first  two  longitudes  at  each  latitude,  the  lat.  and  Ion. 
Indices  and  surface  pressure  are  printed  out,  followed  by  the  12 
sigma  layer  values  of  grid  point  temperature,  winds,  and  specific 
humidity. 


PROGRAM  NAME:  SFPRSl 

ABSTRACT :  This  code  performs  the  univariate  optimum  interpolation  (01) 

analysis  of  terrain  surface  pressure  over  land,  and  a  multivariate 
(with  winds)  01  analysis  of  surface  pressure  over  oceans.  For 
further  discussion  of  the  procedure,  see  Gerlach  (1983). 

AFGL  CYBER  STORAGE: 

SFPRAN/un=NORQUIS  (RECLAIM  from  CC4830) 

AFWL  CRAY  STORAGE; 

Mass  storage:  spx 

On  backup  tape  CJBCK  (AFWL  CYBER):  spx  (fifteenth  file) 

MOS  COMTROL  CARDS  FOR  AFGL  CYBER: 

NORQUIS , T=3000 , CM-27  7  777. 

USER , NORQUIS , NORQUIS . 

CHARGE, 1049, 66 70. 

FTN , SL , R=3 , PL=30000 . 

ATTACH , TAPE1=SP0412 . 

GET,TAPE7=GLAT62. 

ATTACH, TAPE2=F2B0412 . 

GET,TAPE6=TRSST31. 

ATTACH , TAPE4=GR0412 . 

GET ,TAPE5=TRSST62 . 

DEFINE , TAPE9=SFP0412/CT=PU . 

LGO. 

GET,PRASAP. 

ROUTE , PRASAP , DC=IN . 

EOR 

INPUT  FILES: 


TAPEl  Output  (TAPE2)  from  program  PRESP  (see  above). 


TAPB2  Output  from  program  BUFF  (see  Appendix  A;  surface  data  only  are 
used  in  this  code,  Types  3,  7,  8). 

TAPB4  Output  (TAPB4)  from  program  PRBSP  (see  above). 

TAPB5  62  X  62  grid  point  values  of  terrain  surface  height  (a  second 
word  on  each  record  representing  the  corresponding  sea  surface 
temperature  is  not  used) . 

TAPB6  31  X  31  spectral  coefficients  of  terrain  surface  geopotential  (a 
second  record  containing  31  x  31  spectral  coefficients  of  sea 
surface  temperature  is  not  used). 

TAPB7  62  Gaussian  latitudes  (sine  of  latitude)  from  South  Pole  to  North 
Pole  (one  value  per  record). 


ITEMS  REQUIRING  CHANGB  FOR  DIFFBRBNT  TIMBS.  RBSOLUTIONS.  BTC.: 


LINB  NUMBBR(S) 

VARIABLB  OR  CONSTANT 

BXPLAHATION 

25.  26.  43 

2  (OBCALC) 

2,4  (SPTOGPX) 

2,4  (SPTOGPW) 

2  (FFTGP) 

2.3  (LBGSUM) 

2.4  (PMNS) 

2-4(UMVM) 

31  (and  32  when  paired 
with  31) 

Rhomboldal  truncation 
(0-30  waves)  for  spectral 
coefficients 

29,  43 

5  (OBCALC) 

12 

Number  of  sigma  layers  (and 
13  When  it  appears  with  12,  i 
in  array  SI) 

28,  42  (first  ref.) 

4  (OBCALC,  first  ref.) 

8  (FLAGS,  first  ref.) 

32  (SFPRS2.  first  ref.) 

62 

Number  of  evenly-spaced 
longitudes  in  analysis  grid 

28,  42  (second  ref.) 

4  (OBCALC,  second  ref.) 
8  (FLAGS,  second  ref.) 
32  (SFPRS2,  second  ref. 

62  " 

) 

Number  of  Gaussian  latitudes 
in  analysis  grid 
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62-65 

IM0N,IDATE 

Month  (0-12),  date  (0-31), 

IYBAR,IT1MB 

year  (last  two  digits),  and 
hour  (0-23)  of  the  analysis 

56 

10 

Maximum  number  of  observations 

31,  32  (SFPRS2) 

actually  used  to  correct  a 
grid  point  value 

2  (CHLSKY) 

12,  13  (RECALC) 

OUTPUT  FILE: 

TAPE9  62  X  62  record  of  corrected  (analyzed)  surface 

pressures. 

WRITTEH  OUTPUT: 

1.  Listing  of  all  observations  not  used  in  the  analysis 
because  of  missing  data  or  wrong  time. 

2.  Diagnostic  prints  of  values  of  various  parameters  from 
flagging  box  ICB  (see  line  number  20  in  FLAGS). 

3.  Listing  of  box  number,  latitude,  longitude,  JDSI,  p, 
u,  V  residuals,  quality  controls  for  p  and  winds,  and 
model  terrain  height  for  all  observations  used  in  the 
analysis. 

4.  Listing  of  latitude  index,  latitude,  longitude  index, 
longitude,  forecast  pressure  value,  calculated 
correction,  normalized  analysis  error,  and  number  of 
observations  used  in  the  analysis  of  each  grid  point. 
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PROGRAM  MAME;  ADJVBR 


ABSTRACT:  This  program  inputs  the  analyzed  surface  pressure  and  the  grid 

point  forecast  values  of  T,  u,  v,  q  and  interpolates  (linearly 
in  In  p)  these  values  to  the  new  (analyzed)  pressure  values  of 
the  sigma  layers.  It  then  converts  each  interpolated  62  x  62 
field  to  31  X  31  spectral  coefficients,  and  then  re-evaluates 
these  back  on  the  62  x  62  grid. 

AFGL  CYBER  STORAGE: 

PRASAP/un=NORQUIS  (RECLAIM  from  CC4830) 

AFWL  GRAY  STORAGE: 


Mass  storage:  preasap 

On  backup  tape  CJBCK  (AFWL  CYBER):  preasap  (thirteenth  file) 

MOS  COHTROL  CARDS  FOR  AFGL  CYBER: 

HORQUIS , T=1000 , CM=2  7  7  7  7  7 . 

USER , NORQUIS , NORQUIS . 

CHARGE, 1049, 66 70. 

FTH , SL , R=3 , PL=50000 . 

ATTACH , TAPE1=GR0412 . 

ATTACH , TAPE2=SFP0412 . 

GET , TAPE4=GLTWT31 . 

DEFINE , TAPE3»SPN0412/CT-PU . 

DEFINE ,TAPE6=GRN0412/CT»PU . 

GET,BIN=BINCFFT. 

LOAD, BIN. 

LGO. 

GET,ASPRAN. 

ROUTE , ASPRAN , DC=IN . 


IMPUT  FILES: 


TAPEl  Output  (TAPE2)  from  program  PRBSP  (see  above). 

TAPE2  Output  (TAPE9)  from  program  SFPRSl  (see  above). 

TAPE4  31  Gaussian  latitudes  (sine  of  latitude)  from  equator  to  pole,  to 
20  decimal  places,  follotied  by  31  Gaussian  weights  in  same 
orientation  and  accuracy. 


ITEMS  REQUIRIMG  CHAMGE  FOR  DIFFEREMT  TIMES.  RBSOLUTIOMS.  ETC.; 


LIME  WUMBER(S)  VARIABLE  OR  COMSTAMT 


BXPLAMATIOM 


10,  11  (first  ref.), 
14  (first  ref.), 

50 

2,3(FFT1) 

8(SPTOGPG) 

2,4(REORD) 

7(SPT0GP) 


62 


11  (second  ref.),  62 

14  (second  ref.), 

28,  51 

12-14,19,48,49  31  (and  32  when 

2,3(LEGSUN,GQL,LMN)  paired  with  31) 

2(FFT1);  2,4(PMMS) 

2 (PRINT);  3-5,10,15,16 
(SPTOGPG);  2-4, 

8(SPT0GP) 

2-4,  7(UMVM) 


10,12,22,52 

2(LNPINTP) 

3, 8 (SPTOGPG) 

2-4(REORD) 

2,7(SPTOGP) 


12 


16 


DEL 
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NOVRINT 


Humber  of  evenly-spaced 
longitudes  in  the  analysis  grid 


Number  of  Gaussian 
latitudes  in  the  analysis  grid 


Rhomboldal  truncation  for 
spectral  coefficients 
(0-30  waves) 


Number  of  sigma  layers  (and  13 
when  it  appears  with  12,  as  in 
the  array  SI) 


Sigma  layer  depths  (surface  to 
top) 

Switch  to  allow  (0)  or  turn  off 
(1)  vertical  interpolation  to 
new  sigma  layer  pressures 


OUTPUT  FILES: 


TAPES  Spectral  coefficients  of  T,  >)•  D,  q  and  after  T,  u,  v, 

q  grid  fields  were  vertically  Interpolated  to  new  sigma  layer 
pressures. 

TAPE6  62  X  62  X  12  grid  values  of  T,  u,  v,  q  and  62  x  62  values  for 
p^^^  from  evaluation  of  spectral  values  on  TAPES. 

WRITTEM  OUTPUT; 

1.  Sample  output  from  subroutine  LMPIIITP  for  T,  u,  v,  q  at  two  selected 
grid  points. 

2.  List  of  sines  of  Gaussian  latitudes. 


S. 


For  the  first  two  longitudes  at  each  latitude,  the  lat.  and  Ion. 
Indices  and  surface  pressure  are  printed  out,  followed  by  the  12 
sigma  layer  values  of  grid  point  T,  u,  v,  q. 


PROGRAM  MAME: 


ASAPl 


ABSTRACT :  This  code  calculates  level  Z  and  layer  u,  v,  q  residuals  using 

spectral  coefficients  of  T,  n.  D,  q  on  new  sigma  layer 
pressures  and  upper  air  observations  (Types  1,  2,  4,  6a).  It 
then  performs  a  multivariate  01  analysis  for  Z,  u,  v  on  sigma 
layers  and  a  univariate  01  analysis  for  q  on  sigma  layers. 

AFGL  CYBER  STORAGE; 

ASPRAN/un=NORQUIS  (RECLAIM  from  CC4830) 

AFWL  CRAY  STORAGE: 


Mass  storage:  asapoi 

On  backup  tape  CJBCK  (AFWL  CYBER):  asapoi  (seventeenth  file) 

HOS  CONTROL  CARDS  FOR  AFGL  CYBER: 


NORQUIS , T=6000 , CM=305000 . 
USER . NORQUIS . NORQUIS . 
CHARGE, 1049,6670. 
FTN,SL,R=3,PL=99999. 

ATTACH , TAPE1=SPN0412 . 
GET,TAPE7=GLAT62. 

ATTACH , TAPE2=F2B04 12 . 
GET,TAPE5=TRSST31 . 

GET , TAPE10=ESPRER . 

ATTACH, TAPE4=GRN0412 . 

DEFINE , TAPE6  =HMM04 1 2 /CT=PU . 
DEFINE , TAPEl 1 =CE04 1 2/ CT=PU . 
NAP, ON. 

LGO. 

GET,PSTASAP. 

ROUTE , PSTASAP , DC= IN . 

EOR 
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IIIPUT  FILES: 


TAPBl  Output  (TAPB3)  from  program  ADJVBR  (see  above). 

TAPE2  Output  from  program  BUFF  (see  Appendix  A;  upper-air  data  only 
are  used  In  this  code,  Types  1,  2,  4,  6a). 

TAPB4  Output  (TAPB6)  from  program  AOJVBS  (see  above). 

TAPES  31  X  31  spectral  coefficients  of  terrain  surface  geopotential 
(a  second  record  containing  31  x  31  spectral  coefficients  of 
sea  surface  temperature  is  not  used). 

TAPE?  62  Gaussian  latitudes  (sine  of  latitude)  from  South  Pole  to 
North  Pole  (one  value  per  record). 

TAPEIO  Estimated  prediction  error  values  for  Z,  u,  v,  q  on  12 
mandatory  levels  (100-50  mb).  Records  are  structured: 
Z(l,12);  u(l,12):  v(l,12);  q(l,12);  one  record  for  each 
Gaussian  latitude. 


ITEMS  REQUIRING  CHANGE  FOR  DIFFERENT  TIMES.  RESOLUTIONS.  BTC.: 


LINE  NUMBBR(S)  VARIABLE  OR  CONSTANT  EXPLANATION 


21,22,42,125,126,129,130  31  (and  32  when 

4(MASTORl,MASTOR2)  paired  with  32) 

4,5(MAST0R6) 

3,4(FG);  2.4(SPT0GPX, 

SPTOGPW);  2(FFTGP) 

2,3(LEGSUM);  2,4(PMNS) 

2-4(UMVM) 


Rhomboidal  truncation  for 
spectral  coefficients 
(0-30  waves) 


41  (first  ref.) 


62  Number  of  evenly-spaced 

longitudes  in  the  analysis  grid 


41  (second  ref.) 


62  Number  of  Gaussian  latitudes  in 

the  analysis  grid 
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.N 


.1.  JVW  J-  J. 


21,22,24,25,28,29,33,34  12 

40,62,80;  2,3;  8-10 

(MASTORl);  2, 3, 6-8 

(MAST0R4);  16,18,19 

(QCALC);  2,4,5(0BTMZL); 

2,3,4,9-ll(MASTOR2); 

2,3,5,ll-13(HASTOR6); 

2-4,6(CALCRES);3,4,6,7, 

11,12(FLAGS);  2,3,5,6(FG); 

2-4(PTOSIG);  20(TLAZLE); 

23-26,28,30,  31 ,33(ASAP2) ; 

2,4,7,11  (QCOR);  2-5(L0WTIlP) ; 
IKESOBER) 


25 (ASAPl);  3 

,5(L0WTMP) 

25 

Equal  to  2*ML-t-l,  where  ML 
the  number  of  sigma  layers 

28, 109 (ASAP; 

10 (MASTORl); 

85 

Equal  to  7*ML+1 

8(HASTOR4);  19(QCALC); 

5(OBTMZL);  ll(MASTOR2); 

13 (MASTOR6 ) ;  6 (CALCRBS) ; 

12(FLAGS);  6(FG);  33(ASAP2) 

23,37 ,107 (ASAPl) ;  8  (also,  96)  NSHDS,  the  maximum  number  of 

23-28, 40(ASAP2) ;  complete  soundings  that  are 

2-6 (QCOR)  considered  for  contributing  to 

the  correction  of  a  particular 
grid  point  (96sML*lISlIDS) 

10  (also,  20)  NLRHS,  the  maximum  number  of 
individual  height  or  wind 
observations  actually  used  to 
correct  a  grid  point  value  (if 
all  observations  were  winds, 
there  would  be  20,  10  u's  and 
10  v's) 

114-117  IM0H,rDATB,  Month  (0-12,  date  (0-31), 

iyR,ITIMB  year  (last  two  digits),  and 
hour  (0-23)  of  the  analysis 

OUTPUT  FILES: 

TAPE6  62  X  62  X  12  corrected  (analyzed)  values  of  T,  u,  v,  q  in  the 
following  form: 

DO  140  LAT=1,62 
DO  140  LOir=l,62 

140  WRITE(6)(T(K),K=1,12),  (U(K) ,K=1 , 12) ,  (V(K) ,K=1 , 12) , 

(Q(K),K=1,12) 


108 (ASAPl);  27,28 
(first  ref.  only), 
29,30(ASAP2); 
2,6,7(QCOR);  2(CHLSKY); 
12, 13 (RECALC) 


Number  of  sigma  layers  (and  13 
when  it  appears  with  12,  as  in 
array  SIO 


64 


TAPBll  Calculated  corrections  Z  ,  U  ,  V  ,  Q  and  estimated 

c  c  c  c 

analysis  errors,  Z^,  U^,  Vg,  Qg  on  62  x  62  x  12  grid, 
output  In  the  following  manner: 

DO  140  LAT»1,62 
DO  140  L01l=l,62 

140  WSITE(ll)  (Z  (K),K=1,12),  (U  (K),K=1,12), 
c  c 

(V^(k),K=l,12),  (Q^(K),K=1,12),  (Zg(K),K=l,12), 
(Ug(K),K=l,12),  (Vg(K),K=l,12),  (Qg(K) ,K=1 , 12) 

WRITTEM  OUTPUT: 

1.  Number  of  observations  (NUN)  and  timings  (SCND)  from  each  of  the 
MASTORB  subroutines  that  calculate  the  residuals. 

2.  Diagnostic  prints  of  values  of  various  parameters  from  flagging  box 
ICB  (see  line  number  43  in  FLAGS). 

3.  Listing  of  latitude  Index,  latitude,  longitude  Index,  longitude, 
forecast  value,  calculated  correction,  normalized  analysis  error, 
and  number  of  observations  used  in  the  analysis  of  each  variable  of 
grid  points  LATsl,61,  LON=l. 

4.  Listing  of  box  number,  latitude,  longitude,  JDSI,  and  Z,  u,  v,  q 
residuals  and  quality  controls  for  every  fiftieth  observation. 

5.  Listing  of  all  pressure  levels  and  layers,  and  all  Z,  u,  v,  q 
residuals  and  quality  controls  that  were  collected  around  grid 
points  LAT=38  and  LON=55,62. 


PROGRAM  NAME:  CONSPB 

ABSTRACT :  This  program  takes  the  62  x  62  grid  values  of  T,  u,  v,  q  at  12 

sigma  layers  and  which  were  corrected  by  means  of  the  01 

analysis  and  expands  them  Into  spectral  components  for  the  AFGL 
global  spectral  model. 

AFGL  CYBER  STORAGE: 

PSTASAP/un=NORQUIS  (RECLAIM  from  CC4830) 

AFWL  CRAY  STORAGE; 

Mass  storage:  postasapoi 

On  backup  tape  CJBCK  (AFWL  CYBER):  postasapoi  (sixteenth  file) 

NOS  CONTROL  CARDS  FOR  AFGL  CYBER; 

NORQU I S , T= 1 000 , CM=  26  0000 . 

USER . NORQUIS . NORQUIS . 

CHARGE, 1049,6670. 

FTN.SL,R=3. 

ATTACH , TAPE1=MMM0412 . 

ATTACH , TAPE2=SPN0412 . 

GET , TAPE4 =GLTWT3 1 . 

GET, TAPES =TRSST31. 

DEFINE  ,TAPE6=ABIL0412/CT=PU . 

GET,BIN=BIMCFFT. 

LOAD, BIN. 


IMPUT  FILES: 


TAPEl  Output  (TAPK6)  from  program  ASAPl  (see  above). 

TAPE2  Output  (TAPB6)  from  program  ADJVER  (see  above). 

TAPE4  31  Gaussian  latitudes  (sine  of  latitude)  from  equator  to  pole, 
to  20  decimal  places,  followed  by  31  Gaussian  weights  in  same 
orientation  and  accuracy. 

ITEMS  REQUIRING  CHANGE  FOR  DIFFERENT  TIMES.  RESOLUTIONS.  ETC.: 


LINE  NUMBER(S) 

VARIABLE  OR  CONSTANT 

EXPLANATION 

8,9(first  index),  31; 
2,3(FFT1) 

62 

Number  of  evenly-spaced 
longitudes  in  the  analysis  grid 

9  (second  index),  21,62. 

62 

Number  of  Gaussian  latitudes  in 
the  analysis  grid 

10-12,  29-30;  2, 3 (LEG- 
SUM, GQL,LMN)  ,  2(FFT1), 
2,4(PMNS),  2(PRINT, 
FGSHIFT) 

31  (and  32  when 
paired  with  32) 

Rhomboidal  truncation  for 
spectral  coefficients  (0-30 
waves) 

10,12,22,33;  2,3(FGSHIFT) 
2(RMSC) 

12 

Number  of  sigma  layers 

14,  4(FFT1) 

263 

Dimension  of  work  array  WSAVE 
(must  be  >  4*MP+15) 

5 (FGSHIFT) 

DEL 

Sigma  thickness  of  layers 

OUTPUT  FILES; 


TAPE6  31  X  31  spectral  coefficients  of  analyzed  n<  D,  T,  q  on  12 
sigma  layers,  analyzed  arid  in  the  form  of  the 

output  from  program  ASTOGS  (see  Appendix  B). 


WRITTEN  OUTPUT: 

1.  Root  mean  square  values  computed  from  spectral  coefficients  of 
analyzed  n,  D,  T,  q,  and  Values  computed  for  n,  D,  T,  q 

are  for  all  layers,  and  for  a  weighted  average  of  all  layers. 
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END 

FILMED 

8-85 

DTIC 


